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1. INTRODUCTION AND SUMMARY 


Current design philosophy for hydrogen-fueled, scramjet-powered hypersonic 
aircraft results in configurations with strong couplings between the engine plume 
and vehicle aerodynamics. As seen in the sketch below, the aft portion of the 
vehicle will be utilized as the expansion nozzle for the scramjet, and the resulting 
exhaust gas impingement will play a major role in determining the aerodynamic 
forces and mements on the vehicle. In order to assess the impingement effects 
during wind tunnel experiments, the correct thermodynamics as well as aero- 
dynamics of the exhaust gases must be reproduced. 




The method of using the detonation tube simulator to achieve a near exact simulation 
of the exhaust plume to validate cool substitute gases with properties that would allow 
them to behave in a manner similar to the real reacting engine exhaust is discussed 
in Ref. 1. This concept was shown to be valid for the Mach 8 flight case in Ref. 2. 
Present research has extended the proof of concept to the Mach 6 flight case, and has 
obtained data to provide a "bench mark" for comparison of theoretical methods cur- 
rently under development (Ref. 3). 

This report describes the Mach 6 experiments and associated analysis. Pres- 
sure and heat transfer measurements were obtained over a simulated cowl, and over 
an afterbody inclined at both 0 and 20° to the exit nozzle centerline. Scramjet com- 
bustor exit plane conditions were duplicated for the Mach 6 case, and results were 
compared with those obtained for a substitute gas mixture of 50% Argon and 50% 

Freon 13B1. Shock waves were purposely introduced into the exit flow for some runs 
by placing wedges at the cowl exit and nozzle sidewall. Time exposure glow photo- 
graphs were obtained of the shock wave locations. 

Results were excellent. For the shock-free flows the normalized pressure dis- 
tributions obtained in the combustion flow are in almost exact agreement with those 
obtained in the substitute gas flow and with two-dimensional method of characteristics 
calculations where such comparisons were valid. In the cases for which strong shock 
waves were introduced into the flows the agreement among combustion gas, substitute 
gas, and calculations was very good, although the shock wave from the cowl in the sub- 
stitute gas flow appears to impinge on the afterbody slightly upstream of the combus- 
tion gas impingement point. Heat transfer measurements indicate a turbulent boundary 
layer exists over the afterbody with very high heating rates in the region of the nozzle 
exit plane. In our previous tests at M^ = 8 (Ref. 2) heat transfer measurements indi- 
cated the boundary layer was laminar. 

To study further some of the phenomena associated with the Mach 6 flight condi- 
tion, an extension was granted to the original contract for "Validation of Scramjet Ex- 
haust Simulation Technique at Mach 6". The extension had two objectives. The first 
objective was to examine theoretically and experimentally the effects of different sub- 
stitute gas mixtures on the impingement location of a shock wave propagating across 
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the exhaust flow field. The calculated impingement location for various substitute 
mixtures showed very little effect of gas composition. The constancy of shock location 
was verified by experiment with one substitute gas composed of Freon 12 and Argon. 

The other objective was to examine, by experiment, the influence on afterbody 
heating rate distributions of simulated combustion nozzle separators. The experi- 
ments were made using both the substitute and combustion gases for various separator 
positions. The results of these experiments show a wake structure which increases 
the heating rate by approximately 40% over the local heating rate without the separator. 
This result is consistent for both gases but is only evident for X greater than . 78, 
since for more upstream locations the wake is too narrow relative to the spatial reso- 
lution of the instrumentation. 

Use of commercial products or names of manufacturers in this report does not 
constitute official endorsement of such products or manufacturers, either expressed 
or implied, by the National Aeronautics and Space Administration. 

2. DETONATION TUBE SIMULATOR 

The Grumman detonation tube simulator used for the measurement of exhaust 
flows is capable of giving nearly correct chemistry and total enthalpy for the hydro- 
gen/air combustion system throughout the hypersonic cruise flight regime (Ref. 1). 
Reynolds numbers can easily be maintained at flight values by running at elevated pres- 
sures with both combustion and substitute gases. A schematic of the detonation tube 
facility is shown in Fig. 1. It consists of a 6.1 meter long, 7. 62 cm inside diameter 
driver section initially separated by a metallic diaphragm from a 10 meter long, 12.7 
cm inside diameter driven tube. The driven tube is terminated by a supersonic nozzle 
designed to produce a flow that will match the expected exit plane conditions of the 
proposed scramjet engine. The nozzle exhausts into a 1.83 meter diameter, 3.66 
meter long test section. The test section can be evacuated prior to a run to any de- 
sired pressure down to 5 x 10” 5 torr*. The facility has been designed to handle com- 
bustible hydrogen mixtures safely. Various interlocks and leak detection devices are 
incorporated into the automatic gas handling system. Details of this system may be 
found in Ref. 4. A photograph of the facility is shown in Fig. 2. 

Mounted within the test section and mated to the supersonic nozzle is the model 
afterbody. The present experiments were carried out on two different configurations, 

* 1 torr. = 133. 3 Pa 
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Fig. 1 Detonation Tube Simulator Schematic 





Fig. 2 Photograph of Detonation Tube Simulator 



a 20° afterbody surface and a 0° afterbody surface. For each configuration 1/ 20 
scale pressure and heat transfer models were constructed. Details of the models are 
given in Figs. 3 through 12. The backward-running detonation technique was used ex- 
clusively for the current experiments. This technique was developed specifically for 
the simulation of hydrogen/oxygen rocket engine plumes and is described in detail in 
Refs. 1 and 4. The running conditions for the present tests are given in Table 4. 

For running the inert substitute gases in the detonation tube we run the facility 
as a conventional reflected shock tunnel, as described in Ref. 2. The stagnation con- 
ditions of the substitute gas are determined by the strength of the incident shock wave 
and the initial pressure in the driven tube, both of which may be easily controlled. 

The stagnation conditions are given in Table 5 herein. 

DETONATION TUBE/COMBUSTOR NOZZLE DESIGN 

To duplicate the exit flow of the scramjet combustor in the detonation tube facil- 
ity, a two-dimensional supersonic nozzle was employed. The Mach number of the 
supersonic combustion ramjet (SCRJ) exhaust gas at the combustor exit plane for the 
Moo= 6 flight condition is quite low (=»1. 7, see Table 2, Ref. 1). This means the 
throat-to-exit plane area ratios of the detonation tube combustor nozzle must be small 
(^1.4). Because the throat area of the detonation tube nozzle cannot exceed approx- 
imately 10% of cross-sectional area of the detonation tube, the nozzle and consequently 
the scale of the afterbody model involved also must be about l/20 scale, in contrast to 
the 1/8 scale experiments at Moo = 8 of Ref. 2. Our nozzle was designed to duplicate 
the nozzle and model being constructed for substitute gas experiments at NASA Langley 
Research Center so that the two test programs would be directly comparable. Nozzle 
coordinates were supplied to us by Langley personnel and are listed in Fig. 3 together 
with a drawing of the nozzle contour. Since this nozzle was designed for a substitute 
gas mixture of 50% Freon 13B1 and 50% Argon, we checked its suitability for use with 
combustion gas. We calculated the flow field through the nozzle using the M-O-C com- 
puter program of Ref. 5 with detonation tube combustion gas properties generated by 
the computer program of Ref. 6. The difference in the nozzle exit plane Mach number 
distribution between the substitute gas and combustion was insignificant. 
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Fig. 3 Two-Dimensional Contoured Nozzle for Simulating Moo = 6 SCRJ Exhaust Flow 
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Fig. 6 Drawing of 20° Afterbody Pressure Plate 
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NOTE: DIMENSIONS IN CENTIMETERS 


Fig. 7 Drawing of 0° Afterbody Pressure Plate 



POLISHED SURFACE 



NOTE: DIMENSIONS IN CENTIMETERS 

Fig. 8 Drawing of Reflection Plane Pressure Plate 
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Fig. 11 Drawing of Reflection Plane Heat Transfer Plate 
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COWL AND AFTERBODY MODEL DESIGN 


Two afterbody configurations were employed for the present tests: one parallel 
to the nozzle centerline (called the 0° afterbody) and the other inclined at an angle of 
20° to the nozzle centerline. Each configuration consisted of two flat plates: one for 
mounting pressure transducers and one for heat transfer gauges. One cowl configura- 
tion was used with both the afterbody configurations, and was machined as an integral 
part of the nozzle. Provision was made for mounting shock wave generators on the 
cowl lip. With the 0° afterbody a shock generator inclined 3° to the nozzle axis was 
used. With the 20° afterbody the shock generator was inclined 19° to the nozzle axis. 

One side of the model assembly consisted of an instrumented reflection plane, 
extending the entire length of the afterbody from the nozzle exit plane. Instrumenta- 
tion included both heat transfer and pressure transducers. The other side of the model 
consisted of a short side plate extending from the nozzle exit plane to the cowl lip. 

The side plate had provision for shock wave generators to be mounted on its trailing 
edge. The side shock generator was a 15° wedge for use with the 0° afterbody, and a 
20° wedge for use with the 20° afterbody. Figure 4 is a drawing showing the signifi- 
cant dimensions of the nozzle and afterbody model. Figure 5 is a photograph showing 
the model assembled with the 20° afterbody. Figures 6 through 11 are detailed draw- 
ings of the afterbody and reflection plane plates showing all the heat transfer and pres- 
sure transducer locations. 

For the first part of the contract extension, simulated combustor nozzle sepa- 
rators were installed in our nozzle. They extended from the subsonic portion of the 
nozzle to the combustor exit plane, and had a zero-camber, transonic airfoil cross- 
section. They were made in two halves and could be located at any position across 
the width of the nozzle. The nominal positions were at 1/3 and 2/3 of the nozzle width, 
thus dividing the 3:1 (width to height) combustor exit plane into three 1:1 segments. 
Figure 12 is a sketch of these separators and shows their installation in the nozzle. 
Figure 13 is a photograph of the installation. 

To study the near- field effects of the wakes of these separators, additional in- 
strumentation had to be added to the 20° afterbody heat transfer model. The new in- 
strumentation consisted of 13 additional heat transfer gauges and 7 new pressure 
transducers. Their locations are shown shaded in Fig. 9. 
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DISTANCE ALONG PLATE SURFACE (INCHES) 
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PRESSURE 
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Fig. 12 Sketch Showing Nozzle Separators and Their Installation 
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One side effect of the separators was that their finite thickness (0. 533 cm; 

0.210 inch) decreased the throat area of our combustor nozzle and consequently in- 
creased the nozzle exit Mach number. All heat transfer data taken with sectors in 
place had to be corrected for these changes in nozzle exit conditions. The correction 
procedure is described in the Experimental Results Section (5). 

The second part of the contract extension consisted of measuring the cowl shock 
wave impingement location on the 20° afterbody in different substitute gas flows. In 
order to ensure that the shock wave and its impingement location were completely 
within a two-dimensional flow region, the nozzle side plate, which originally extended 
to the cowl lip (X/Y = 3.12), was extended beyond the expected impingement location 

O 

to X/Y 3 = 6. 62. 

INSTRUM ENT ATION 

Model heating rates were measured with thin film transducers, Ref. (7). These 

gauges, when used without taking into account changes in substrate properties with 

2 

temperature, have a range of 0 to 50 cal/cm sec and an accuracy of ±7% in this range. 
When corrections are made for changes in substrate properties (Ref. 8), the useful 
range is increased to 500 cal/cm sec with an accuracy of ±15%. For this work all 
but a few of the combustion heating rates measured on the 20° afterbody and about 
half of the points on the 0 afterbody were in the 500 cal/cm sec range. Those heat- 
ing rates above 500 cal/ cm sec were corrected by extrapolating the correction curve 
and are therefore subject to greater error . All data taken using the substitute gas are 
in the most accurate range. Scatter of the data for the substitute gas runs and most of 
the combustion runs should be assigned to turbulence and not gauge accuracy. 

Model pressures were measured with semiconductor strain gauge pressure 
transducers while tunnel operating pressures and impact pressure were obtained with 
piezoelectric pressure transducers. Both gauges are described in Ref. 2 and their 
performance characteristics are outlined in Table 1. The spatial resolution for de- 
tecting shock location is superior for the heat transfer gauge. This advantage plus the 
lower cost per unit were the reasons for using heat gauges for determining shock lo- 
cations before committing pressure instrumentation. 
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TABLE 1 MEASUREMENT ACCURACY 


Tunnel Parameter 

Accuracy 

Range 

P^ Tunnel Driven Initial Pressure 

P, Tunnel Driver Initial Pressure 

4 

Pj_ Stagnation Pressure 

U. Incident Wave Speed 

me 

Prp^ Impact Pressure 

Q Heat Transfer Rate 

±689 N/m 2 
±1.38 x 10 5 N/m 2 
±3.10 x 10 5 N/m 2 
±30 m/sec 
±6.89 x 10 4 N/m 2 
±7% 

0 -►6. 89 x 10 5 N/m 2 

0-+-1.38 x 10 7 N/m 2 

0-^1. 38 x 10 7 N/m 2 

0-*-900 m/sec 

0-*-6.89 x 10 6 N/m 2 
2 

0-*-50 gm cal/cm sec 


All model data were recorded on a 35-channel FM multiplex recorder system 
after suitable signal conditioning. Data were digitized one channel at a time by an 
A/D converter and a time average of the signal during the tunnel test time was ob- 
tained from this digital information. A more detailed description of the data acqui- 
sition and reduction system appears in Ref. 2. Major components of the system are 
outlined in Table 2. 

SPATIAL RESOLUTION OF HEAT TRANSFER GAUGE IN THE Z DIRECTION 

Our heat transfer gauges have an active area defined by a span of 0. 762 cm 
(perpendicular to the flow) and a chord of 0. 159 cm (parallel to the flow). Each mea- 
sured heating rate was the average over this area. Therefore, strong gradients in the 
spanwise z direction produce measured heating rates smaller than the peak value. For 
the experiments performed under the contract extension it was not possible to reori- 
ent the instrumentation (rotate the gauges 90°); with the existing gradients in the y 
direction such reorientation would have been of doubtful value. In order to assess the 
probable impact of limited gauge resolution, the following discussion may be helpful. 

Given a rectangular perturbation in heating rate such as: 


- q max 

[-*-% GAUGE WIDTH 


GAUGE WIDTH 
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TABLE 2 INSTRUMENTATION AND DATA RECORDING EQUIPMENT 



Pressure Transducers 


Manufacturer 

Model No. 

Max Range 

Kuhlite 

VQS-250-200 

1.38 x 10 6 N/m 2 

Kuhlite 

VQL- 250-20 

1.38 x 10 5 N/m 2 

Kuhlite 

VQH-250-5 

3.45 x 10 4 N/m 2 

Kistler 

603 L 

1. 03 x 10 7 N/m 2 

Kistler 



607 

4.13 x 10 8 N/m 2 

Data Recording and Reduction 

Honeywell 

7600 tape recorder 

1. 5 MHz bandwidth 

D.C.S. 

Gov-3 VCO's 

±2v, ±32 T kHz 

D. C.S. 

GMA-5, line drivers 

- 

D. C.S. 

GFD-14, discriminators 

±10v, ±16 kHz 

Hewlett- Packard 

A/D converter 

12 bit, 50>usec/word 

Hewlett- Packard 

21 MX computer 

- 

Texas Instruments 

700 A SR terminal and 
cassette tape drive 
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a linear averaging of heating rate over the gauge gives a result 50% lower than the 
peak. Obviously, this average is in greater error for narrower perturbations. For 
high heating rates as in the combustion gas case the problem is further complicated 
by the temperature dependence of the gauge substrate properties, which makes the 
equations relating the heating rate to time and temperature nonlinear. Averaging of 
such high rates gives an even lower result for the peak heating rate. However, this 
additional error is small compared to the error caused by spatial resolution alone. 

FLOW VISUALIZATION 

Glow photographs, i.e. , those taken by the light of the flow, gave low contrast 
prints (Figs. 14 and 15) when no shock generators were used. When the top shock 
generator or top and side generators were used a bright region in the position of the 
shock could be seen (Figs. 16 and 17). These photos, taken with Polaroid PN type 55 
film with a 135 mm lens or with a 254 mm lens, show shock locations for the combus- 
tion runs, but the contrast is still quite low and the visible plume appears to fill the 
entire test section. Schlieren photos could not be made in this configuration since the 
reflection plane is in the schlieren optical path. 


3. DETONATION TUBE RUNNING CONDITIONS 

SHOCK TUBE MIXTURE DEVELOPMENT 

Any one of a number of initial fuel mixtures may be used to achieve a proper 
simulation of the scramjet exhaust. For each of these mixtures there is an appro- 
priate initial pressure and incident shock speed that will yield the correct pressure, 
enthalpy, and chemistry in the stagnation region. 

Our first attempts used an initial fuel mixture of hydrogen and air. Runs made 
using this mixture had a large time delay between the arrival of the incident shock at 
the end wall and the formation of a detonation wave. There was burning of the mixture 
at shock arrival, seen by the large increase in pressure and noise, but that the detona- 
tion wave was not established until much later (see Fig. 25). Efforts to shorten this 
delay by increasing incident shock velocity had no consistent effect until the incident 
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Fig. 16 Glow Photograph, Side View, 20° Afterbody Model with Top (Cowl) ShockWave Generator 






wave became strong enough to cause ignition of the incident shock, what we call a 
front-running detonation. We also tried starting the detonation with two spark plugs 
rather than one, with the hope that a more uniform ignition source would shorten the 
delay time. This also had no consistent effect. At this point we decided to investigate 
other mixtures. 

In previous scramjet simulation of the = 8 condition (Ref. 2), we used a fuel 
which had all of the oxidizer in the form of nitrous oxide. For the = 6 condition 
this fuel required high initial pressures and low incident shock velocities which make 
shock tube operation difficult. Nevertheless, a series of runs was made with nitrous 
oxide substituted for various amounts of oxygen in the fuel. None of these runs pro- 
duced less delay time than hydrogen and air. Since having the oxidizer in a different 
form failed, we tried supplying some of the required hydrogen contained in ammonia. 

We started with a mixture that had all of the hydrogen in the form of ammonia 
and found it produced a front-running detonation when subjected to a shock of the proper 
velocity for correct simulation. Since the ammonia sensitized the reaction, a series 
of tests were made to determine what proportion of ammonia and hydrogen was needed 
to reduce the ignition delay without igniting the incident wave. This line of testing es- 
tablished the fuel mixtures described previously and was used for all rims to obtain 
data for the combustion case. 

COMBUSTION GAS (0 = 1.0) 

A complete description of the operation, performance, and computational pro- 
cedures involved in determining the detonation tube simulator running conditions is 
given in Ref. 1. All the combustion gas data presented herein were obtained using the 
backward-running detonation technique. 

As shown in Fig. 7 of Ref. 1 there are any number of gas mixtures that are 
apparently eligible to be used for simulation. However, we encountered considerable 
difficulty in getting many of these mixtures to ignite properly. Some would ignite very 
early, during the passage of the incident wave, and others would ignite very late, after 
the passage of the reflected wave. Repeatability was also a problem, with things like 
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the condition of the shock tube wall appearing to be a factor affecting the rate of igni- 
tion (e.g. Ref. 9). Section 5 of this report gives the details of these experimental 
ignition problems. 

We eventually found one mixture which behaved repeatably and for which ignition 
delay was tolerable (Z.1. 0 msec. ). This mixture was composed of 

0. 220 H 2 + 0. 157 0 2 + 0. 533 N £ + 0. 063 NH g + 0. 007 Ar 

To determine the state of the test gas after the delayed ignition we carefully measured 
the incident and reflected shock wave speeds and formulated a one-dimensional wave 
diagram consistent with the wave speed measurements, the pressure time- history 
measurements, and previous models of delayed reflected ignitions (e.g. Ref. 10). 
Figure 18 shows this wave diagram. 

Figure 19 is an oscilloscope picture showing the detonation tube stagnation pres- 
sure time history for this running condition. (Note the difference in the pressure time- 
history between this running condition with a short ignition delay and that with zero de- 
lay as shown in Fig. 13a of Ref. 2. ) We were then able to calculate the state of the gas 
in regions 6 and 7 using the computer program of Ref. 6, by assuming that a Chapman- 
Jouget detonation started in region 5 and then propagated into region 2. Table 3 sum- 
marizes the results of these calculations and gives the state of the gas in the various 
regions involved. Note that the region 6 gas is about 60 K cooler than ideal, while the 
region 7 gas is about 160 K hotter than ideal. This is illustrated in the Mollier dia- 
gram sketch shown in Fig. 20. In practice these two regions would tend to mix, due 
to diffusion, turbulence, etc. , rather than remain as two discrete regions, and the re- 
sulting flow out of the nozzle should be extremely close to the desired, ideal testing 
condition. 

The measured static pressure at the exit of the nozzle, P„, and on the flat por- 

6 9 

tion of the cowl (before the 6 degree bend) was 1.730 x 10 N/m (251. psia). An os- 
cilloscope picture showing this measurement is presented in Fig. 21. Based on the 
measured Pp./P 3 ratio, the ideal total enthalpy, and assuming an isentropic, equili- 
brium expansion, we computed (Ref. 6) the exit plane Mach number to be 1.706. Table 
4 presents the significant running conditions obtained by the procedure outlined above. 
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Fig. 18 Distance-Time Wave Diagram Illustrating Delayed Start of Reflected Detonation Wave 



Fig. 19 Oscilloscope Record of Detonation Tube Stagnation Pressure with Delayed Reflected Ignition 
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Fig. 20 Mollier Diagram Sketch of Detonation Tube Stagnation Conditions 
1 atm = 101.3 kPa 
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Fig. 21 Oscilloscope Record of Combustor Nozzle Exit Plane Static Pressure 
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TABLE 3 EQUILIBRIUM THERMODYNAMIC STATE PROPERTIES IN 
BACKWARD- RUNNING DETONATION TUBE FOR SIMULATING 1/20 
SCALE SCRJ EXHAUST AT Moo = 6 (0 = 1.0) 

Region 1 0. 21967 H 2 + 0. 15691 0 2 + 0. 55348 N 2 + 0. 06276 NH^ + 0. 00719 AR 

P, = 1.364 atm, T 1 = 298. K, U = 996.34 m/sec 

1 IS 

Region 2 same gas 

P_ = 10. 043 atm, T = 633. K, U 0 = 708. 84 m/sec 
Region 5 same gas 

P,, = 43. 657 atm, T c = 1001. K, U„ = 404. 88 m/sec 

5 D ±A 

Region 6 stagnant, reacted test gas 

0. 63410 N 2 + 0. 31727 H 0 + 0. 01867 H 2 + 0. 00978 0H + 0. 00783 Ar + 

0. 00568 + 0. 00420 0^ + 

P„ = 80.3 atm, T c = 2918. K, H_ = 237.7 cal/gm, S a = 2.3671 cal/gm K 

b b b b 

Region 7 stagnant, reacted test gas 

0. 62586 N 9 + 0. 29744 H 0 + 0. 02995 H + 0. 01709 0H + 0. 00775 Ar + 

ut Ci 

0. 00916 N0 + 0. 00646 0 £ + 

P^ = 80.3 atm, T ? = 3140. K, H^ = 390. cal/gm, S^ = 2.4169 cal/gm K 


TABLE 4 COMBUSTION GAS RUNNING CONDITIONS 
<£= 1.0, or+ 0 = 4°, 1/20 Scale 



Ideal 
(Ref. 1) 

Based on 
Measured P5/P3 

P 5 (N/m 2 ) 

8. 086 x 10 6 

8. 136 x 10 6 

H 5 (cal/gm) 

278.0 

278. 0 

t 5 (K) 

2982. 0 

2982.0 

P 3 (N/m 2 ) 

1.861 x 10 6 

1. 730 x 10 6 

T 3< K > 

2375. 0 

2343. 0 

m 3 

1. 657 

1.706 

73 

1.2122 

1.2144 

Rey. No. 3 (m -1 ) 

5.02 x 10 7 

4.895 x 10 7 
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SUBSTITUTE GAS 


Thermodynamic properties of the substitute gas mixture used in the present tests were 
computed by the method detailed in Appendix A of Ref. 2 and are listed in Table 5. The 
particular mixture chosen for these tests was 50% Argon and 50% Freon 13B1 with a stag- 
nation temperature of 533. 3 K. This mixture was the gas planned to be used in the con- 
current Langley Test Program (see page 4) and for which the two-dimensional contoured 

nozzle was designed. To achieve full scale Reynolds number at the combustor exit plane a 

5 2 

stagnation pressure of 5. 516 x 10 N/m (80. psia) was required. 

The detonation tube facility was operated in the reflected shock tube or shock tunnel 
mode (Ref. 1) in order to generate the required stagnation conditions for the substitute gas 
runs. 

The incident shock wave strength required to produce the stagnation temperature was 
computed from ideal gas shock equations using an average value of 7 across the incident and 
reflected shock waves. Figure 22 is an oscilloscope picture showing the reflected shock 
stagnation pressure in the shock tube. Note' that as reported in our previous work at the 
M oo = 8 flight condition (Ref. 2), the same reflected shock bifurcation phenomena occurs; 
this results in a somewhat gradual rise to the final stagnation pressure level. A detailed 
explanation of this phenomena is given in Ref. 2. 

For these substitute gas runs the static pressure at the exit of the nozzle and on the 

5 2 

flat portion of the cowl, P^, measured 1.434 x 10 N/m (20. 8 psia). From the measured 
Pg/Pg ratio and the calculation procedure described earlier, we determined the Mach num- 
ber at the exit plane of the nozzle to be 1. 63. This agrees very closely with the value calcu- 
lated for our nozzle using the two-dimensional method-of-characteristics program (Ref. 5). 

One of the objectives of the contract extension was to study the sensitivity of the cowl 
shock wave impingement location with variations in substitute gas properties. To do this we 
selected five additional substitute gas mixtures and/or stagnation temperatures and calculated 
the nozzle and afterbody flow fields with each of these gasses. The additional gas mixture/ 
temperature combinations selected were 

50% Argon + 50% Freon 13B1 @ T5 = 477. 8 K 

60% Argon + 40% Freon 13B1 @ T5 = 533. 3 K 

60% Argon + 40% Freon 13B1@ T5 = 477. 8 K 
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60% Argon + 40% Freon 12 @ T5 = 533. 3 K 
60% Argon + 40% Freon 12 @ T5 = 477. 8 K 


The results of these calculations are discussed in the next section. Tables 6 through 10 
give the thermodynamic properties of these mixtures , determined by the procedure detailed in 
Appendix A of Ref 2. 

i 

TABLE 5 THERMODYNAMIC PROPERTIES OF 50% ARGON - 

50% FREON 13B1 SUBSTITUTE GAS MIXTURE ATT = 533. 3°K 

o 


Mach Number 

y 

T(K) 

cal 

Cp( K m K> 

P/Po 

5.4882 

1.3695 

111.1 

0.07799 

3. 6027E-C4 

5.0921 

1.3451 

127.8 

0.08022 

6. 1281E-04 

4.5523 

1.3127 

155.6 

0.08834 

1. 3580E-03 

4. 1114 

1. 2877 

183.3 

0. 09419 

2. 7686E-03 

3.8795 

1.2753 

200.0 

0.09749 

4. 1148E-03 

3.4039 

1. 2519 

238.9 

0.10459 

9. 6545E-03 

3. 0470 

1. 2364 

272. 2 

0. 11006 

1. 8794E-02 

2.8268 

1.2278 

294.4 

0.11342 

2. 8510E-02 

2. 6166 

1. 2203 

316. 7 

0.11656 

4. 2430E-02 

2.4138 

1.2137 

338.9 

0.11949 

6. 2070E-02 

2. 2155 

1. 2080 

361. 1 

0. 12223 

8. 9396E-02 

2.0194 

1. 2028 

383.3 

0. 12478 

0. 12693 

1.8720 

1. 1994 

400.0 

0.12658 

0. 16366 

1. 6724 

1. 1952 

422. 2 

0.12883 

0. 22723 

1.4117 

1. 1906 

450.0 

0.13143 

0. 33697 

1. 2443 

1. 1882 

466. 7 

0. 13287 

0.42343 

1. 0623 

1. 1859 

483.3 

0. 13424 

0.52907 

0. 9969 

1. 1852 

488.9 

0. 13468 

0. 56915 

0. 8555 

1. 1838 

500.0 

0. 13553 

0. 65749 

0.4852 

1. 1813 

522.2 

0. 13713 

0.87146 

0. 0000 

1. 1801 

533.3 

0. 13789 

1.0000 
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TABLE 6 THERMODYNAMIC PROPERTIES OF 50% ARGON - 

50% FREON 13B1 SUBSTITUTE GAS MIXTURE ATT = 477. 8 K 

o 


Mach Number 

y 

T(K) 

~ , cal . 

P/Po 

5. 0403 

1.3695 

111.1 

0.07799 

7.3299E-04 

4. 5515 

1.3378 

133.3 

0.08333 

1. 4737E-03 

4. 0574 

1.3072 

161. 1 

0.08955 

3.2047E-03 

3.6462 

1. 2833 

188.9 

0.09531 

6. 4432E-03 

3.2893 

1. 2644 

216. 7 

0. 1006 

1. 2202E-02 

2. 9692 

1. 2490 

244.4 

0.1055 

2. 2032E-02 

2. 6743 

1. 2364 

272.2 

0. 1101 

3. 8238E-02 

2.3959 

1. 2258 

300.0 

0. 1142 

6. 4176E-02 

2. 1802 

1.2186 

322. 2 

0. 1173 

9. 5084E-02 

2. 0203 

1. 2137 

338. 9 

0. 1195 

0. 1 2628 

1.8602 

1.2093 

355.6 

0.1216 

0. 16626 

1. 7525 

1.2066 

366. 7 

0. 1229 

0.19879 

1. 6981 

1. 2053 

372.2 

0. 1235 

0.21709 

1. 6433 

1. 2041 

377.8 

0. 1242 

0.23687 

1. 5319 

1.2017 

388. 9 

0. 1254 

0.28129 

1.4172 

1. 1994 

400.0 

0.1268 

0.33298 

1.2980 

1. 1972 

411.1 

0. 1277 

0.39293 

1. 1066 

1. 1942 

427.8 

0.1294 

0. 50090 

0.9661 

1. 1924 

438. 9 

0. 1304 

0. 58681 

0. 8084 

1. 1906 

450.0 

0.1314 

0.68558 

0. 6201. 

1. 1890 

461. 1 

0.1324 

0. 79885 

0. 3547 

1. 1874 

472. 2 

0. 1333 

0.92845 

0. 0000 

1. 1866 

477.8 

0. 1338 

1. 00000 
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TABLE 7 THERMODYNAMIC PROPERTIES OF 60% ARGON- 
40% FREON 13B1 SUBSTITUTE GAS MIXTURE AT T q = 533. 3 K 


■ Mach Number 

7 

T(K) 

cal 

P/Po 

5.1299 

1.4003 

111.1 

0.083216 

7. 9151E-04 

4.6384 

1.3716 

133.3 

0.087802 

1. 5238E-03 

4.3317 

1.3534 

150.0 

0. 091092 

2. 3726E-03 

4.0626 

1.3376 

167.7 

0. 094257 

3. 5765E-03 

3.7473 

1.3193 

188.9 

0.098282 

5. 9346E-03 

3.5357 

1.3074 

205.6 

0.10116 

8. 4590E-03 

3.2786 

1.2936 

227.8 

0.10481 

1. 3192E-02 

3.0436 

1.2816 

250.0 

0.10825 

2. 0014E-02 

2.7732 

1.2689 

277.8 

0.11226 

3. 2614E-02 

2.5215 

1.2581 

305. 6 

0.11596 

5. 1520E-02 

2.2361 

1.2472 

338. 9 

0.12003 

8. 6106E-02 

2.0062 

1.2395 

366. 7 

0.12311 

0.12879 

1.9152 

1.2367 

377.8 

0.12427 

0.15042 

1.8242 

1.2341 

388.9 

0.12538 

0.17513 

1. 7329 

1.2317 

400. 0 

0.12646 

0.20330 

1.6870 

1.2305 

405.6 

0.12698 

0.21880 

1.5945 

1.2282 

416. 7 

0.12800 

0.25291 

1.5006 

1.2261 

427.8 

0.12899 

0.29154 

1.4048 

1.2241 

438. 9 

0.12994 

0. 33521 

1.3062 

1.2222 

450.0 

0.13085 

0. 38444 

1.1511 

1.2195 

466. 7 

0.13216 

0.47005 

1.0407 

1.2178 

477.8 

0.13299 

0.53592 

0.9827 

1.2170 

483.3 

0.13340 

0.57176 

0.7190 

1.2140 

505.6 

0.13495 

0.73677 

0.4487 

1.2119 

522.2 

0.13604 

0.88628 

0.0000 

1.2106 

533.3 

0.13674 

1.00000 
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TABLE 8 THERMODYNAMIC PROPERTIES OF 60% ARGON - 
40% FREON 13B1 SUBSTITUTE GAS MIXTURE AT T q = 477.8 K 


Mach Number 

y 

T(K) 

^ , cal , 
^gm 

P/Po 

4. 6249 

1.4389 

111.1 

0.07799 

1. 6860E-03 

4.0885 

1.3911 

138. 9 

0. 08461 

3. 6120E-03 

3.5824 

1.3492 

172.2 

0.09191 

8. 0182E-03 

3. 1020 

1.3138 

211.1 

0. 09960 

1. 8188E-02 

2.8035 

1. 2944 

238. 9 

0. 10459 

3. 0914E-02 

2. 5290 

1. 2785 

266. 7 

0. 10919 

5. 0676E-02 

2.3212 

1. 2678 

288. 9 

0. 11260 

7. 3595E-02 

2. 1201 

1. 2585 

311. 1 

0.11580 

0.10504 

1. 9227 

1.2504 

333.3 

0.11878 

0.14760 

1. 7750 

1.2450 

350.0 

0. 12089 

0.18873 

1. 7255 

1. 2433 

355.6 

0. 12156 

0.20450 

1.6759 

1. 2416 

361. 1 

0. 12223 

0.22141 

1. 5756 

1. 2385 

372. 2 

0. 12353 

0.25894 

1.4217 

1. 2341 

388. 9 

0. 12539 

0.32563 

1.2050 

1. 2289 

411. 1 

0.12773 

0.43766 

1. 0886 

1. 2264 

422. 2 

0. 12883 

0. 50535 

1.0275 

1.2253 

427. 8 

0. 12937 

0. 54250 

0. 9639 

1. 2242 

433.3 

0. 12990 

0.58203 

0. 7508 

1.2210 

450.0 

0. 13143 

0.71610 

0.5760 

1. 2190 

461. 1 

0.13240 

0.81982 

0.3295 

1. 2171 

472. 2 

0. 13336 

0. 93644 

0.0000 

1.2162 

477. 7 

0. 13379 

1. 00000 
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TABLE 9 THERMODYNAMIC PROPERTIES OF 60% ARGON - 
40% FREON 12 SUBSTITUTE GAS MIXTURE AT T q = 533. 3 K 


Mach Number 

y 

T(K) 

c p& 

p/p o 

5.2088 

1.3973 

111. 1 

0.09661 

6. 5371E-04 

4. 6104 

1.3581 

138. 9 

0. 10419 

1. 4765E-03 

4. 1386 

1.3281 

166.7 

0.11119 

3.0162E-03 

3. 6745 

1.3006 

200. 0 

0.11886 

6. 4694E-03 

3. 1040 

1. 2707 

250.0 

0.12897 

1. 7694E-02 

2.7751 

1.2558 

283.3 

0. 13484 

3. 2271E-02 

2.5210 

1. 2458 

311. 1 

0.13925 

5. 1458E-02 

2. 2788 

1.2373 

338. 9 

0. 14325 

7. 9879E-02 

2. 0438 

1. 2301 

366. 7 

0.14685 

0. 12109 

1.8577 

1.2251 

388. 9 

0. 14948 

0.16632 

1. 7644 

1. 2229 

400.0 

0. 15072 

0.19399 

1.7175 

1.2218 

405.6 

0. 15131 

0. 20927 

1. 6704 

1.2208 

411. 1 

0. 15190 

0.22559 

1. 5753 

1. 2188 

422.2 

0. 15303 

0.26158 

1.4786 

1. 2169 

433.3 

0.15412 

0.30246 

1.3795 

1. 2151 

444.4 

0. 15516 

0.3488 

1. 1707 

1. 2119 

466. 7 

0. 15711 

0.46026 

1. 0582 

1.2104 

477.8 

0. 15802 

0.52677 

0. 9991 

1.2097 

483.3 

0. 15847 

0.56305 

0. 8042 

1.2077 

500. 0 

0. 15974 

0.68520 

0. 5609 

1. 2058 

516. 7 

0.16094 

0. 82973 

0. 3209 

1. 2047 

527. 8 

0. 16170 

0. 94016 

0. 0000 

1.2041 

533.3 

0.16207 

1.00000 
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TABLE 10 THERMODYNAMIC PROPERTIES OF 60% ARGON - 
40% FREON 12 SUBSTITUTE GAS MIXTURE AT T q = 477.8 K 


Mach Number 

y 

T(K) 

_ cal 
P ^m 

P/Po 

4.7917 

1.3973 

111.1 

0.09661 

1. 2410E-03 

4.2217 

1.3581 

138. 9 

0. 10419 

2. 8030E-03 

3.6877 

1.3230 

172.2 

0. 11252 

6. 5438E-03 

3. 2505 

1. 2967 

205. 6 

0. 12006 

1. 3835E-02 

2. 8728 

1.2764 

238. 9 

0. 12686 

2. 7181E-02 

2.6429 

1.2653 

261. 1 

0.13100 

4. 1262E-02 

2.4251 

1.2558 

283. 3 

0. 13484 

6. 1262E-02 

2.2158 

1.2476 

305.6 

0.13840 

8. 9183E-02 

2.0114 

1.2405 

327. 8 

0. 14170 

0.12756 

1. 9101 

1. 2373 

338. 9 

0. 14325 

0.15164 

1.8087 

1. 2343 

350.0 

0. 14473 

0.17958 

1. 7069 

1. 2314 

361. 1 

0. 14616 

0.21189 

1. 6040 

1.2288 

372.2 

0. 14753 

0.24916 

1.4463 

1.2251 

388. 9 

0.14948 

0.31573 

1.3376 

1. 2229 

400. 0 

0. 15072 

0.36827 

1.2247 

1.2208 

411.1 

0. 15190 

0.42826 

1. 1060 

1.2188 ! 

422.2 

0.15303 

0.49658 

0.9788 

1. 2169 : 

433.3 

0. 15412 

0. 57419 

0.8389 

1. 2151 

444.4 

0.15516 

0.62614 

0. 5843 

1. 2127 

461. 1 

0. 15664 

0. 86599 

0.3341 

1. 2112 

472.2 

0. 15757 

0.93502 

0.0000 

1.2104 

477.8 

0. 15802 

1.0000 
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4. AFTERBODY FLOW FIELD PREDICTIONS 


For the calculations presented all j combustion gas thermodynamic data were de- 
rived from the computer program of Ref. 6. Substitute gas thermodynamic data were 
generated by the procedure described in Appendix A of Ref. 2. The flow fields, in- 
cluding the presence of shock waves, were calculated using the two-dimensional, MOC 
program of Ref. 5. The geometry of the afterbody and cowl surfaces is fiven below. 


The coordinate system is shown in Figs. 

23 and 24. 


Cowl 


Y 

= 1.0 

0 < X < 1.11 

Y 

= 0. 4204 x 2 - 0. 933 X + 1. 518 

1.11 < X < 1.235 

Y 

= 0.1051 X + 0.8768 

1.235 < X <3.12 (no wedge) 

Y 

= 0.1051 X + 0.8768 

1.235 < X <2.6 (wedge for 20° afterbody) 

Y 

= -0.3443 X + 2.0453 

2.6 < X < 4. 0 

Y 

= 0.1051 X + 0.8768 

1.235 < X <2.6 (wedge for 0° afterbody) 

Y 

=-0.05241 X + 1.2863 

2.6 < X < 3.5 

20 

0 Afterbody 


Y 

o I 

= (0.1736 - X ) 2 - 0.4167 

0 < X < 0.1425 

Y 

= -0.3640 X + 0.02674 

0.1425 < X < 21.67 

0° 

Afterbody 


Y 

= 0.0 

0 < X <21.67 


where Y = Y/Y 0 , X = X/Y 0 , and Y q is the height of the combustor exit. For our model Y 0 

O u O 

1. 524 cm (0. 6 inch). For the full scale vehicle we assumed Y g = 30.48 cm (1. 0 foot). 

As a result of these calculations we found that the basic afterbody flow field 
and its pressure distribution were relatively insensitive to small changes in combus- 
tion gas thermodynamic properties and combustor exit Mach number. This point was 
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Fig. 23 Two-Dimensional Flow Field for 20° Afterbody with and without Shock Generator. Combustion Gas, a + P 
4°, Moo = 6 
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I 


EXPANSION WAVES 




clearly illustrated in Fig. 15 of our previous contractor report (Ref. 2). The same 
was true of small changes in substitute gas properties. We also found that the pres- 
sure distribution for the substitute gas chosen for the present tests (50% Freon 13B1 + 
50% Argon) and for the combustion gas were in extremely close agreement. There- 
fore, in Figs. 23 and 24 we show only one flow field calculation for each geometry 
involved and in Figs. 26 through 58, where the experimental data are also plotted, 
only one theoretical, two-dimensional pressure distribution is shown. Figure 25 
shows an oscilloscope record of detonation tube stagnation pressure with burning be- 
hind reflected shock and overly delayed detonation. 

Because this model had a reflection plane extending the entire length of one 
side and a side plate extending to the end of the cowl on the other side, almost all of 
the data were obtained in two-dimensional flow regions and no theoretical corrections 
for side expansion effects are shown. 
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Theoretical heat transfer distributions for the afterbody without externally 
generated shocks were calculated for us by the NASA Langley Research Center using 
the computer code of Price and Harris (Ref. 11) for turbulent boundary layers, and 
assuming combustion gas products with the pressure distributions shown in Figs. 26 
through 58. The calculations are shown in Figs. 59 through 73. Calculations 
assuming a laminar boundary layer predicted heating rates an order of magnitude 
lower than the turbulent cases. 

As part of the contract extension, afterbody flow fields with the cowl shock 
impinging on the 20° afterbody were computed for the five additional substitute gases 
listed on pages 13 and 14. The objective here was to predict the sensitivity of the shock 
impingement location with variations in substitute gas properties and then verify the 
results experimentally. Since we would be testing with a fixed contour combustor 
nozzle, variations in substitute gas properties would produce variations in the com- 
bustor nozzle exit plane Mach number and?. These latter variations were taken 
into account when doing the afterbody calculations. The calculated flow fields were 
all very similar to those shown in Fig. 23, and the predicted changes in shock im- 
pingement location were insignificant. Table 11 presents a summary of these 
calculations which shows how little the major parameters of the flow field are 

altered with the different substitute gases. 

Since the calculations predicted insignificant changes in flow parameters in the 
different substitute gases, we selected only one of the five mixtures for the experi- 
mental tests. That gas was the mixture of 60% argon + 40% Freon 12 @ T_ = 477.8 K 

5 

and is referred to as substitute gas no. 2 in the Experimental Results Section (5). 



Fig. 25 Oscilloscope Record of Detonation Tube Stagnation Pressure 
with Burning Behind Reflected Shock and Overly Delayed 
Detonation 
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P/P3 



Fig. 40 Pressure Distribution: 20° Afterbody, with Top 
and Side Shock Generators, Row E 










Fig. 41 Reflection Plane Normalized Pressures (P/P3) 20° Afterbody, No Shock Generators 
(Combustion Data Above Substitute Data) 



Fig. 42 Reflection Plane Normalized Pressures (P/P3) 20° Afterbody, with Top (Cowl) Shock Generator 
(Combustion Data Above Substitute Data) 














Fig. 50 Pressure Distribution: 0 Afterbody, with Top 
(Cowl) Shock Generator, Row B 


A COMBUSTION GAS 
□ SUBSTITUTE GAS ■ 1 
20 MOC PREDICTION 





Fig. 52 Pressure Distribution: 0° Afterbody, with Top 
(Cowl) Shock Generator, Row D 


A COMBUSTION GAS 
□ SUBSTITUTE GAS = t 
?D MOC PREDICTION 






Fig. 56 Pressure Distribution: 0° Afterbody, woth Top 
and Side Shock Generators, Row C 


A COMBUSTION GAS 
□ SUBSTITUTE GAS • t 

JDMOC PREDICTION <DOES' 

NOT INCLUOE EFFECT 
OF SIOE GENERATOR! 



Fig. 58 Pressure Distribution: 0° Afterbody, with Top 
and Side Shock Generators, Row E 
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Fig. 59 Reflection Plane Normalized Pressures (P/Pg) 0 Afterbody, No Shock Generator 
{Combustion Data Above Substitute Data) 


- „ .76' > +40^ 


Fig. 60 Reflection Plane Normalized Pressures (P/P3) 0° Afterbody, with Top (Cowl) Shock Generator 
(Combustion Data Above Substitute Data) 













VO/E ' VE/E 



Fig. 65 Heat Transfer Rate Distribution, 20° Afterbody, 
No Shock Generators, Row D 



Fig. 66 Heat Transfer Rate Distribution, 20° Afterbody, 
No Shock Generators, Row E 



Fig. 67 Heat Transfer Rate Distribution, 20° Afterbody, 
with Top (Cowl) Shock Generator, Row A 



Fig. 68 Heat Transfer Rate Distribution, 20° Afterbody, 
with Top (Cowl) Shock Generator, Row B 







Fig. 69 Heat Transfer Rate Distribution, 20° Afterbody, 
cn with Top (Cowl) Shock Generator, Row C 



Fig. 70 Heat Transfer Rate Distribution, 20° Afterbody, 
with Top (Cowl) Shock Generator, Row D 



Fig. 71 Heat Transfer Rate Distribution, 20° Afterbody, 
with Top (Cowl) Shock Generator, Row E 



Fig. 72 Heat Transfer Rate Distribution, 20° Afterbody, 
with Top and Side Shock Generators, Row A 













A - COMBUSTION GAS 
□ - SUBSTITUTE GAS # 1 

TURBULENT PREDICTION (HARRIS) 

(DOES NOT ACCOUNT FOR SHOCK 


• 


A 




GENERATORS) 

2 DMOC PREDICTED 

COWL SHOCK LOCATION 




Fig. 73 Heat Transfer Rate Distribution, 20 Afterbody, with Top and Side Shock Generators, Row B 





TABLE 11 COMPARISON OF VARIOUS FLOW FIELD PARAMETERS 
FOR DIFFERENT SUBSTITUTE GASES AND COMBUSTION GAS 



1/8 Scale 

0.5F-13B1 + 0.5Ar 

0.4F-12 + 0.6Ar 

0.4F-13B1 +0.6Ar 


Combustion Gas 

+ j3 = 4° 

Substitute Gas 
No. 1 

Substitute Gas 
No. 2 

Substitute Gas 
No. 5 


Desired 

Nozzle 

To (K) 


Equilibrium 

Equilibrium 

477.8 

533.3 

477.8 

533.3 

477.8 

533.3 

m 3 

1.657 

1.702 

1.705 

1.702 

1.713 

1.711 

1.716 

1.713 

m 

0.3792 

0.3730 

0.3773 

0.3812 

0.3667 

0.3708 

0.3620 

0.3670 

1911 

2.220 

2.270 

2.269 

2.259 

2.299 

2.289 

2.309 

2.299 

P R.S. 

P 3 

0.5364 

0.5619 

0.5651 

0.5639 

0.5713 

0.5698 

0.5735 

0.5718 

m r.s. 

1.996 

2.011 

2.014 

2.013 

2.016 

2.016 

2.016 

2.016 

X R.S. 

y 3 

5.6435 

5.6417 

5.6424 

5.6410 

5.6502 

5.6475 

5.6488 

5.6494 


NOTE: Subscript R.S. refers to the reflected shock region 
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5. EXPERIMENTAL RESULTS 


PRESSURE DATA 

All afterbody pressure data are presented by rows and normalized to P^, the 
nozzle exit static pressure, which is different for combustion and substitute gas. 
Figures 26-30 show the pressure distribution for the 1/20 scale 20° afterbody with 
no shock generators. These data show excellent agreement between the substitute 
and combustion gases. There is also good agreement with the 2-D prediction for 
gauge rows close to the reflection plane. The gauge rows farthest from the reflection 
plane deviate most from the 2-D prediction, particularly at large values of X/Y 

O 

(length along afterbody normalized by combustor nozzle exit height). This is the ex- 
pected behavior since it is in this region that deviation from 2-D behavior will be 
greatest. Reflection plane data for this configuration, Fig, 41, are presented in 
schematic form. The combustion data are the upper numbers while the substitute 
gas data are the lower. These data show the same good agreement between substitute 
and combustion gases as the afterbody data. 

Pressure distributions for the 20° afterbody with the cowl generator and both 
test media are shown in Fig. 31-35. The agreement between substitute and combus- 
tion gas data is good but the shock arrival is slightly farther downstream for the 
combustion gas. The substitute gas shock arrival is in close agreement with the 2-D 
prediction and both pressure distributions show the same kind of deviation from 2-D 
behavior as in the no-generator case. The reflection plane data, Fig, 42, show the 
same agreement between test media and indicate the shift in cowl shock location for 
the combustion gas. 

Figures 36-40 show the 20° afterbody pressure distributions with both shock 
generators and both gases. The 2-D distribution on these figures is included only 
for reference, since it does not allow for the side shock. The data for this condition 
appear to show poorer agreement between gases than for the previous configurations, 
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but these differences are due to differences in cowl shock location for different gauge 
stations. The intersection of the two shocks gives rise to the highest pressures 
measured on the 20° afterbody; however the absolute level of these pressures is the 
average pressure on the gauge diaphragm and not a true "point” maximum. Since all 
pressure measurements are averaged over the 3.175 mm diameter diaphragm, flow 
fields with strong pressure gradients, i.e. shocks, will yield measured pressures 
different from the actual pressure at a point. Both these data and the reflection 
plane data. Fig. 43, indicate a merging of the side and cowl shocks at their point of 
intersection since there is no localized pressure increase after the cowl shock in row 
A or on the reflection plane. 

The data in Figs. 44-60 are for the 0° afterbody. These data show reasonably 
good agreement between the two test media and also between both gases and the 2-D 
prediction. There is more scatter for these data than for the 20° afterbody because 
of the very high heating rates for the combustion gas. Even though the pressure trans- 
ducer diaphragms were protected by an ablative silicone coating, many gauges were 
destroyed in these tests, particularly those close to the nozzle exit plane. Data from 
gauges that were destroyed by these heating rates could still be reduced provided the 
gauges were intact for the greater part of the steady flow time (3 msec). Destruction 
of a gauge is easily recognized on the signal. Pressure distributions for the combus- 
tion gas on the 0° afterbody are not as complete as for the 20° afterbody because of 
severe gauge attrition; however Figs. 49-53 indicate the cowl shock for the combus- 
tion gas is farther downstream than the substitute gas cowl shock, as was the case 
with the 20° afterbody. 

Cowl pressure data (Fig. 61) show very good agreement between experiments 
and prediction as well as good agreement between the two gases. Pressure data in 
tabular form as a function of X/Y^ for both the 20° and 0° afterbody models in both 
combustion and substitute gas are presented in Tables 12 through 15. 

During the contract extension pressure measurements with substitute gas no. 2 
were made on the afterbody to determine shock location. No separators were in- 
stalled for these runs. The afterbody pressures appear in Table 16 and show good 
agreement with the No. 1 substitute gas pressures since both the shock location and 
normalized pressure levels agree with the previous tests. Cowl pressures with 
substitute gas no. 2 were measured on the nozzle separator center lines just aft 
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TABLE 12 P/P3 20 DEG AFT BODY COMBUSTION GAS 


ROW A 


NO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 


X--'Y3 

P.--P3 

X.-'Y3 

P--P3 

X-'Y3 

P.--P3 

2 . 36 

0. 37 

2 . 36 

0. 36 

£ . 36 

0. 37 

3 . 92 

0. £9 

3. 92 

0.29 

3.92 

0 ■ 3£ 

5.10 

0 . 2 0 

5. 1 0 

0. 19 

5.10 

0. 24 

5.98 

0 . 1 5 

5 . 98 

0 . 43 

5 . 98 

0.53 

6. 86 

0. 13 

6. 86 

0 . 4 0 

6 . 8 6 

0. 52 

7.74 

0 . 11 

f m • r' 4 

0. 26 

7.74 

0. 38 

8 . 6c‘ 

0. 08 

8 . 6£ 

0.15 

8 . 6£ 

0. 31 

9.5 0 

0 . 1 0 

9 . 5 0 

0.15 

9.50 

0 . £2 



ROW 

E 



NO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 



P.--P3 

X Y3 

P.-P3 

X--Y3 

P-P3 

5. 39 

0. 17 

2. 95 

0. 32 

£ . 95 

0. 36 

6 . £6- 

0 . 1 3 

4.51 

0. £4 

4.51 

0. £5 



5 . 39 

0. 19 

5. 39 

0. 48 



6 . £6 

0 . 4 0 

6 . iz‘6 

0 . 68 



ROW 

c 



NC GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 


X •■■■’ V 3 

P--P3 

X Y3 

P--P3 

X Y 3 

P.--P3 

1.18 

0. 37 

1.18 

0. 35 

1.18 

IJ . 36 

3.92 

0. 36 

3. 92 

0 . 36 

5. 68 

0 . 6 1 

4 . 8 0 

0.21 

4 . 8 0 

0 . £ 0 

6 . 57 

0. 47 

5 . 6 8 

0. 14 

5 . 68 

0. £8 



6 . 57 

0 . 1 0 

6 . 57 

0 . 36 





ROW 

D 



NO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 


X Y 3 

P.--P3 

X Y 3 

P.--P3 

X Y 3 

P.-P3 

1 . 77 

0 . 39 

1 . 77 

n m *-|7 

1 . 77 

0 . 38 

2. 95 

0 . 35i 

2 . 95 

0. 34 

£ . 95 

0 . 33 

5. 0 0 

0. 23 

5. 0 0 

0. ££ 

5 . 0 0 

0. 81 

5.88 

0 . 1 6 

5. 88 

0. 34 

5 . 88 

0.37 

6. 76 

0 . 1 0 

6 . 76 

0. £9 

6.76* 

0. 37 

7. 64 

0. 07 

7.64 

0. 18 





ROW 

E 



NO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 

X--Y3 

P--P3 

X--Y3 

P--P3 

X-‘Y3 

P.-P3 

4. 32 

0.23 

4. 32 

0. £3 

4. 32 

0.47 

6. 08 

0. 08 

6 . 08 

U .26 

6 . 08 

0. 12 

6. 96 

0. 06 

6.96 

0.17 

6 . 96 

0. 09 
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TABLE 13 P/P3 20 DEG AFT BODY SUBSTITUTE GAS 




ROW 

A 


NO GENERATORS 


TOP GENERATOR 

EOTH GENERATORS 

X-'Y3 

P--P3 

X Y 3 

P-P3 

X.--Y3 P.--P3 

£ . 3 6 

0 . 36 

5.10 

0 . £6 

5.10 0.38 

3. 9£ 

0. £9 

5 . 93 

0. 48 

5.93 0.5£ 

5 . 1 0 

0 . £ 0 

6 . 36 

0 . 3 6 

6.86 0.39 

5. 93 

0 . 14 

7. 74 

0. £7 

7.74 0.35 

6 ■ 3 6 

0. 13 

9 . 5 0 

0 . 15 

9 . 5 0 0 . £ 0 

7.74 

0 . 1 0 




9 . 5 0 

0. 09 






ROW 

B 


NO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 

X Y 3 

P--P3 

X s' Y 3 

P.--P3 

X--Y3 P-P3 

£ a 95 

0 . 33 

£.95 

0 . 38 

£.95 0.34 

4.51 

0 . £ 0 

4.51 

0. £3 

4.51 0. £3 

5 . 39 

0 . 1 5 

5 . 39 

0. £6 

5.39 0.45 

6 . £6 

0.11 

6 . £6 

0 . 39 

6. £6 0. 55 



ROW 

c 


NO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 

X Y 3 

P.--P3 

X Y 3 

P--P3 

X ■•••'>■' 3 P.--P3 

3 . 9£ 

0. £9 

4 . 3 0 

0 . £ 1 

4 . 8 0 0. 51 

4 . 3 0 

0 . 1 8 

5 . 68 

0. 34 

5.63 0.53 

5 . 6 8 

0 . 1 3 

6.57 

0. £7 

6.57 0 . £ 7 

6 . 57 

0. 09 






ROW 

B 


NO GENERATORS 


TOP GENERATOR 

EOTH GENERATORS 

V...-Y3 

P--P3 

X •■'Y3 

P--P3 

X.-Y3 P.--P3 

£.95 

0. 35 

£ . 95 

0. 38 

£.95 0.37 

5. 0 0 

0.17 

5. 0 0 

0. 18 

5 . 0 0 0 . 59 

Ej O i“j 

0.15 

r o i - . 
,_i ■ O iZ> 

0. 16 

5.88 0.£5 

6.7b 

0.11 

6 . 76 

0 . 3 1 

6.76 0 . £ 0 



ROW 

E 


NO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 

X • Y 3 

P.--P3 

X V 3 

P.--P3 

X.--Y3 P--P3 

4. 3£ 

0 . £ 0 

4 . 3£ 

0. £4 

4 . 3£ 0.5£ 

6. 03 

0. 09 

6, 08 

0. £6 

6.03 0. 1£ 

6 . yfct 

0. 06 

6 . 96 

0.17 

6 . 96 0. 09 
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TABLE 14 P/P3 0 DEG AFT BODY SUBSTITUTE GAS 





— 

RD l.i.l 

H 


MD GEMERRTDRS 



TDP 

GEMERRTDR 

EDTH GEMERRTDRS 

f V 3 

P.--P3 

V 

Y3 


P--P3 

X7Yy 

£ . 9.8 

0 . 7£ 

iZ. m 

92 


0.71 

£ . 92 

3 . 75 

0 . 5 0 

Z' 

75 


0. 5£ 

3 . 75 

4 . 58 

0 . 5 0 

4. 

58 


1.10 

4 . 58 

5 . 63 

0 . 44 

5. 

63 


0 . 7 0 

5 . 63 

6.67 

0 . £ 1 

6 • 

67 


0. £9 

6 . 67 

7.50 

0 . 14 

f ■ 

5 0 


0 . 1 9 

7.50 

O O |”j 

0 . 12 

8 . 

■~» “i 


0 . 14 





— 

ROW 

E 


MO GEMERRTDRS 



TOP 

GEMERRTDR 

EDTH GEMERRTDRS 

X Y 3 

P--P3 


V 3 


P--P3 

X..-Y3 

3 . 96 

0. 46 

1 . 

£5 


1 . 0 £ 

1.25 

4.79 

0 . 5£ 

“• 

96 


0 . 65 

3 . 96 

5.33 

0 . 37 

4. 

79 


1 . £6 

4. 79 

6 . 67 

0 . £ 0 

5- 

O 


0 , 64 

Cj » 2 j •“« 

7.50 

0 . 13 

b • 

67 


0. £4 

6.67 

9.17 

0. 09 

7 m 

5 0 


0 . £6 

7 . 5 0 



9 . 

17 


0. 09 

9.17 




— 

RDW 

f; 


MD GENERATORS 



TOP 

GEMERRTDR 

EDTH GEMERRTDRS 

! :: !-’V3 

F.--P3 

y 

Y3 


P--P3 

X,-Y3 

c! . i' 1 

0 . 89 

C m 

71 


0 . 78 

£.71 

3 . 75 

0 . 53 

■Z* > 

7 cj 


0. 54 

cr 

Z* . i 

4.53 

0. 58 

4. 

cr 


1.15 

4 . 58 

5 . 42 

0 . 5 0 

c 

4£ 


0 . 86 

5. 4£ 

6 . £5 

0 . £ 1 

& m 

£5 


0 . 37 

6 . £ 5 

7.50 

0 . £ 0 

r** ■ 

5 0 


0. £3 

7 . 5 0 




— 

RDW 

D 


NO GEMERRTDRS 



TDP 

GEMERRTDR 

EDTH GEMERRTDRS 

X ■■■' Y 3 

P--P3 


Y 3 


P--P3 

X.-"Y3 

0 . 83 

1 . 04 

0 . 

83 


1 • 0 £ 

0 . 83 

3. 12 

0 . 63 

m 

1 £ 


0 . 6 1 

3 . 1 £ 



4. 

O t‘ 


1 . £ 0 

4.37 

5. £1 

0 . 59 

5. 

£1 


1.01 

5 . £ 1 




— 

RDW 

E 


MD GEMERRTDRS 



TDP 

GEMERRTDR 

EDTH GEMERRTDRS 

>v-"Y3 

i*i O 

P--P3 

|“i "7 i"t 

v .. 

Y • -! 

l“l C 1 


P.-P3 


3.75 

U ■ i u 

0 . 54 

C. m 

l_l r» 

75 


ij . f y 
0. 54 

c! . 08 
3.75 

4. 79 

0 . 53 

4. 

79 


0 . 93 

4.79 


P.--P3 

0.71 

0 . 31 

1. SO 
0 . 69 
0 . 34 
0.31 
0. 16 


P.-P3 

1 . 04 
1 . 04 
1 . £3 
0. 61 
0 . 33 
0. £4 
0 . 1 0 


P.--P3 

0 . 79 

1.11 

1 . £9 
0. 94 
0. 37 
0 . £ 1 


P--P3 
1 . 01 
0 .59 
1.35 
0 . 93 


p.-'F3 

0. 77 

1 . 16 
0. 75 
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.TABLE 15 P/P3 0 DEG AFT BODY COMBUSTION GAS 


ROW R 

NO GENERATORS TOP GENERATOR BOTH GENERATORS 


X.-'Y3 

P.--P3 

X Y 3 

P.--P3 

‘y* 

P.--P3 

£. OS 

0 . 74 

4 . 53 

1 . 0£ 

4 . 53 

0 . 00 

4. 53 

0. 5£ 

6 . 67 

0 . 3 0 

6.6-7 

0. £8 

6 . 6- 7 

0. £4 

7.50 

0 . £ 0 

7.. 5 0 

0. £4 

7 . 5 0 

0. £4 







ROW 

B 



NO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 


X Y o 

P.--P3 

X--'Y3 

P.-P3 

X.-Y3 

P.--P3 

1 . £5 

1 . 0 0 

3 . 96 

0. 53 

3 . 9 6 

0. 73 

3 . 96 

0. 50 

4. 79 

1 . 13 

4 . 79 

0. 9£ 

4. 79 

0.59 

5. 33 

0. 55 

cr -i 

i . 

0 . 49 

6 . 67 

0. £3 

6 . 67 

0. 34 

6 . 67 

0 . 36 

7.50 

0. 14 

7 . 5 0 

0. 10 

? . 5 0 

0 . £ 0 



ROW 

f: 



NO GENERATORS 


TOP GENERATOR 

BDTH GENERATORS 


X Y - J 

P--P3 

X ,-Y 3 

P--P3 

X ' Y3 

P.--P3 

£.71 

0 . 8 0 

4.53 

1 . 0£ 

3 , 75 

0. 34 

3.75 

0 . 59 

5. 4£ 

0. 39 

5 . 4£ 

0. 76 

4.5S 

0 . 59 

6 . £5 

0. 4£ 

6 . £5 

0. 37 

5. 4£ 

0. 43 

7 . 5 0 

0. 18 

7.50 

0. 16 

6. £5 

U ■ c! r 





7.50 

0. 16 






ROM D 

NO GENERATORS TOP GENERATOR BOTH GENERATORS 

X-- Y3 P.-P3 X--Y3 P.-'P3 X 3 P-"P3 

5.31 0.57 7.5 0 0.18 3.13 0.57 


NO GENERATORS 

X.--Y3 P.--F'3 

£.09 0.60 

3.75 0.57 


KCI.I E 

TOP GENERATOR BOTH GENERATORS 

X-"Y3 P--'P3 X’-'t'S P P 

7.50 0.13 3.75 1.30 

4.73 0. 3£‘ 



TABLE 16 P/P3 20 DEG AFT BODY SUBSTITUTE GAS 2 

NO GENERATOR TOP GENERATOR 



R.-'P 3 

P--P3 

- - 33S 

. £31 


-.459 

. 153 

. 143 

1 . 3£0 

. 097 

. 3 1 9 

-1 . 3£8 


. 1 04 

-.495 


. 094 
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of the nozzle exit and, when scaled by P5 and compared to the scaled pressure with- 
out the separators, did not show any significant changes. Note that the active gauge 
area is large compared to the wake for this X/Yg location. 

Cowl pressure measurements were attempted using the combustion gas with 
nozzle separators. However, the first run in this configuration had a very large 
ignition delay and all of the gauges were destroyed without yielding useful data. These 
were the only gauges provided in this pressure range and we could not continue with 
this measurement. 

IMPACT PRESSURE MEASUREMENTS 

We were unable to obtain meaningful impact pressure measurements because of 
the large size of our probe relative to this nozzle. We suggest future free-stream 
measurements of this type be made with a Laser Doppler Velocimeter system, since 
it does not require a probe which disturbs the flow. 

HEAT TRANSFER DATA 

Distributions of heating rates for both substitute gas and combustion products on 
the 20° afterbody show fairly good agreement with the turbulent predictions calculated 
for us by NASA using the program of Price and Harris (Ref. 11). All runs made with 
combustion products burned out many heat transfer gauges because the high heating 
rates gave rise to gauge temperatures high enough to melt the solder holding the leads 
to the gauge. Gauge destruction was particularly bad close to the nozzle exit for the 
0° afterbody case, and in areas where the cowl and or side shock crossed the model. 
Even in cases where a gauge was destroyed by the run, useful measurements were 
obtained as long as the gauge was intact for the greater part of the test time. The 
high attrition rate curtailed our testing of the 0° afterbody with combustion gas; how- 
ever, enough points were taken on the 0° afterbody to show the cowl shock location 
and to supply an increased basis for correlation of heating rates between combustion 
and substitute gases. 

Ail of our testing was performed at a single Reynolds number, matched to the 
flight of a vehicle with a nozzle height of 30. 5 cm (1 ft), for Mach 6 flight at a nomi- 
nal altitude of 24,207 m (79,400 ft). Our tests were l/20th scale. Because the 
Reynolds numbers are matched and the gas properties and Mach number distributions 
are matched for the combustion case, one can expect measured Stanton number (i.e. , 
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q/pu A T) to be the same as flight if the ratio of wall temperature to flow total 
temperature is matched (Ref. 1). Wall temperature effects will be discussed below. 
For matched Stanton number (in the high Mach number, high gas temperature cases 
of interest) one calculates flight absolute heating rates by dividing measured model 
heating rates by the ratio of flight to model dimensions to correct for the density 
ratio. This is because Reynolds number invariance requires the product p D to be 
invariant between model and flight. 

In order to show a comparison between substitute and combustion gas data, the 

levels shown in Figs. 74 and 62-78 are levels measured on the model which have been 

2 

normalized by 388 cal/cm sec for the combustion data and the Price and Harris pre- 

2 

diction. The substitute gas data have been normalized by 2. 71 cal/cm sec. The 

2 2 
388 cal/cm sec value was chosen as a convenient scale factor while the 2. 71 cal/cm 

2 

sec value was computed from the 388 cal/cm sec divided by an average of combustion 
data points, multiplied by an average of substitute data for the same locations. 

Figures 62-66 show the distribution for the 20° afterbody with no generator. 

The agreement between the combustion data and prediction is quite good considering 
the expected accuracy of heat transfer data in this range and the nature of turbulent 
heating predictions in general. The substitute gas data agree with both the combus- 
tion data and the prediction in shape (the average levels are forced to agree by the 
scaling used). Heat transfer distributions for the 20° afterbody with the cowl genera- 
tor are given in Figs. 67-71. In this case the solid curve is intended for reference 
only since it does not consider any shocks. The arrival of the cowl shock for the 2-D 
prediction (X = 5. 63) is very close to the substitute gas arrival (X = 5. 68) and some- 
what ahead of the combustion gas arrival (X = 5. 87). In Figs. 72-75 the distributions 
for the 20° afterbody indicate a "merging” of the cowl and side shock at their inter- 
section. The intersection of these shocks (Fig. 73) gives rise to the highest com- 
bustion gas heating rate measured on the 20 afterbody model (1358 cal/cm sec). 

The heat transfer data for the 20° afterbody are presented in tabular form as a func- 
tion X/Y in Tables 17 and 18. 

O 

The heating rates for the 0° afterbody obtained in both substitute and combustion 
gases appear in Tables 19 and 20. Data for the 0° afterbody with the combustion gas 
are somewhat sparse due to the aforementioned gauge destruction problems. However, 
the cowl shock impingement location for the combustion case is again seen to be aft 
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Fig. 74 Heat Transfer Rate Distribution, 20° Afterbody, with 
Top and Side Shock Generators, Row C 



Fig. 75 Heat Transfer Rate Distribution, 20° Afterbody, with 
Top and Side Shock Generators, Row D 
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Fig. 76 Reflection Plane Normalized Heating Rates (Q/QA), 20 Afterbody, No Shock Generators 
(Combustion Data Above Substitute Data) 
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Fig. 77 Reflection Plane Normalized Heating Rates (Q/QA), 20° Afterbody, with Top (Cowl) Shock Generator 
(Combustion Data Above Substitute Data) 
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\TABLE 17 Q/QA 20 DEG AFT BODY COMBUSTION GAS 
'{QA = 388. cal/cm^ sec) 




ROW 

A 



HD GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 


X,Y3 

Q.--QA 

X Y 3 

Q,- QA 

X--Y3 

Q--QH 

■t; Cj --i 

0. 76 

cr •-! 

0 .61 

5 . 1 0 

0. 50 

5 . 1 0 

0 . 48 

5 . 1 0 

0. 45 

5 . 88 

0.71 

5.38 

0 . 35 

5 . 1 0 

0.51 

7.45 

1. 19 

7.45 

0 . 3 0 

cr 

._i , or, 

0. 95 

9 . 0 1 

0. 70 

9. 01 

0 . £ 1 

7.45 

0. 59 





9 . 0 1 

0. 39 





ROW 

B 



MO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 


X-"Y3 

Q.--QA 

X Y 3 

q.--'Qh 

X--Y3 

Q.-'QA 


0. 88 

4. 51 

0. 65 

4.51 

0 . 53 

5. 89 

0. 45 

5 . £9 

U • 6 6 

5 . £9 

0. 54 

6. 08 

0. 88 

5 . £9 

0 . 57 

6 . 08 

3 . 5 0 

7.45 

0. 84 

6. 08 

1 . 04 

7.45 

0. 89 



7. 45 

0 . 99 





ROW 

C 



NO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 


X V 8 

Q.--QA 

X ■■■'' Y 3 

Q..- QR 

X Y 3 

Q.-'QA 

3.11 

0 . 88 

4.71 

0 . 5 0 

4.71 

0.41 

3. 98 

0 . 8 1 

5 . 49 

0. 37 

5.49 

1 . £ 0 

4.71 

0. 53 

6 . c! 7 

1 . £9 

b ■ c r 

1 . 36 

5. 49 

0 . 39 

7. 45 

0. 78 



7.45 

U . £ 1 







ROW 

I, 



NO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 


X-'Y3 

Q.--QA 

X-"Y3 

Q.-'QA 

X Y 3 

Q--QA 

3. 34 

0. 79 

4 . 9 0 

0 . 49 

4 . 9 0 

1 . 15 

4. 18 

0. 88 

5. 69 

0 . 4 0 

5.69 

0 . 78 

5 . 69 

0. 39 

5 . 69 

0 . 38 

6 . 4? 

0. 93 



6.47 

1.11 

7.45 

0. 55 



ROW 

E 



NO GENERATORS 


TOP GENERATOR 

BOTH GENERATORS 


X ■■■' Y 3 

Q.-'QA 

X Y 8 

Q.-'QH 

X Y 3 

Q-'QA 

cr 

0 . 86 

5 . 1 0 

0 . 56 



5. 88 

0. 44 

5 . 88 

0. 48 





6 . 6 6 

0 . 53 
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TABLE 18 Q/QA 20 DEG AFT BODY SUBSTITUTE GAS 
(QA = 2, 71 cal/cxn 2 sec) 


NO GENERATORS 


'm' y 

Y3 

Q.-'QA 

o 

53 

0 . 68 

c* 

.j a 

10 

0 . 48 

5. 

y y 

0.45 

6 » 

67 

0 . 85 

9 m 

01 

0 . 8 0 


NO GENERATORS 



, : v- i 

Q.- 

QA 

8. 95 


ft 9 

■ ■ f o 

0. 

76 

4.51 


48 

5.89 


48 

6. 08 


j CT 

7.45 

0 . 

c 


NO GENERATORS 



X.--Y3 

Q 

QA 

3.11 

0. 

94 

3 . 9 8 

0 . 

68 

4.71 


59 

5.49 

l"l 

44 

6.87 


4 0 

7.45 


c r 


NO GENERATORS 



X--Y3 

Q..- 

UA 

3. 34 


76 

4.18 

0 . 

55 

4 . 9 0 



5 . 69 

0. 

Z'C. 

6. 47 


31 

7. 45 

o. 

a -z* 


NO 

GENERATORS 

vy 

Q.-QH 


CO 
..i . 

0 . 9 0 

Cj 

*"i o 

0.51 


ROW 

H 

TOP GENERATOR 

X Y 3 

Q.--QA 

3. 53 

0. 65 

5. 1 0 

0. 38 

5 . 88 

0 . 76 

7.45 

0, 59 

9. 01 

0.33 

ROW 

E 

TOP GENERATOR 

X Y 3 

Q.-QH 

8. 95 

0 . 7 0 

•"I* "7 ■' Z ■ 

0 . 66 

4.51 

0. 49 

5. 89 

0. 38 

6. 08 

0.97 

7.45 

0 . 6 7 

ROW 

f: 

TOP GENERATOR 

X •••'' 3 

Q--QR 

3.11 

0. 95 

3.98 

0 ■ 64 

4. 71 

0. 54 

5.49 

0 . 5 0 

6 ■ c! r* 

n ■ 77 

7.45 

0. 58 

ROW 

D 

TOP GENERATOR 

X--V3 

Q.-'QA 

3 . 34 

0. 78 

4.18 

0. 48 

4 . 9 0 

0. 38 

5.69 

0 . 6 1 

6.47 

0. 68 

7.45 

0.59 

ROW 

E 

TOP GENERATOR 

X--Y3 

Q.-'QA 

6 . 67 

0. 43 

7. 45 

0. 59 


BOTH GENERATORS 


X Y 3 

Q.-QA 

3 . 53 

ij . 66 

5.10 

0. 48 

5 . 88 

0 . 96 

7.45 

1. 08 

9 . 0 1 

0 . 65 

BOTH GENERATORS 


y / y ; 

Q--QA 

8 .95 

0. 81 

O . i‘" o 

0 . 85 

4.51 

0 . 55 

5.89 

0. 84 

6 . 08 

1 . 48 

7.45 

0. 95 

BOTH GENERATORS 


X--Y3 

Q.-'QA 

3 . 1 1 

0. 98 

3 . 98 

0 . 6 0 

4.71 

0. 95 

5.49 

1 . 86 

6.87 

1 . 15 

7.45 

0. 59 

BOTH GENERATOR': 


X--Y3 

Q.-'QA 

3 . 3 4 

0. 88 

4. 18 

0 . 77 

4 . 9 0 

1.57 

5 . 69 

1. 17 

6 . 47 

0 . 77 

7 . 45 

0. 57 

BOTH GENERATORS 


X.--'Y3 

Q.-'QA 

3 . 53 

1.87 

5.88 

0 . 75 
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TABLE 19 Q 0 DEG AFT BODY COMBUSTION GAS 


ROW E 


GENERATORS 



TDP 

GENERA TDP 

X V 3 

CfiL 

X V 3 

CRL 

C- - d 9 

6 03. 


4. 79 

537. 00 

3 . 1 3 

595. 


5. 63 

693. 00 

3 . 96 

•jij r* a 

fi f | 

6 . 46 

374. 00 

4.79 

539. 


7 . 3 9 

341 . 00 

5. 63 

347. 

0 0 



6 . 46 

363 . 

0 0 



6 ■ 46 

374 . 

0 0 



7 . 39 

133. 

0 0 



7 . 3 9 

181. 

0 0 






— 

RCIW C 

GENERATORS 



TOP 

GENEERTHR 

X ■■■' V 3 

CHL 

X Y 

CRL 

3.50 

6 1 3 . 

|-| M 

5. 00 

367. 00 

5 . 0 0 

494 . 

|-| i”r 

nr i-i 

. j B o . "i 

619. 00 

nr ,-j •-> 

._» a |“| .Jj 

336 . 

j"i i"i 

r’ « r’ 1 

1 8 0 . 0 0 

7.71 

134. 

u j-j 



7.71 

133. 

0 0 






— 

RDM n 

GENERATORS 



THP 

GENERATOR 

X 't'3 

CPL 

X V >! 

CRL 

1.67 

651 . 


4.37 

431 . 00 

3.71 

603. 

M 

5. 31 

96 0. 00 

3.54 

534. 

0 0 

6. 35 

337. 00 

4.37 

330. 




5.31 

517. 




6 . 35 

356 . 

0 0 



6 . 35 

366. 




7.50 

139. 




7 . 5 0 

139. 

n n 
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TABLE 20 

Q 0 DEG AFT BODY SUBSTITUTE GAS 



ROW 

fi 


MD GENERATORS 


TDP GENERATOR 

EDTH GEMEPRTDRS 

X V 3 

CRL 

, : V Y 8 

CRL 

X--Y3 CRL 



£ . 03 

y - yy 

£.03 4.76 



i — ■ 7* u. 

3.43 

£ . 9£ 4.35 



ROW 

B 


MD GENERATORS 


TDP GEMERRTDR 

EDTH GEMEPRTDRS 

X V • •! 

CRL 

X, 3 

CRL 

X Y 3 CHL 

£. £9 

4. 39 

£ . £ 9 

4. 39 

£ . £ 9 4.73 

3. 13 

4. 13 

3 . 1 3 

3 . 3 9 

1 . 1 o — 1 . r r 

3. 96 

£ . 3 0 

3.96 

£. 93 

• -! . 9 6 £ . 96 

4 . 79 

3. £6 

4. 79 

6.44 

4.79 4.39 

5 . 6 3 

£ . 39 

5.63 

3 . 4 0 

5.63 3.67 

6 » 46- 

1 . 63 

6.46 

£ . 3 6. 

6.46 £ . U7 

P. £9 

1 . 01 

7 . £9 

1 .£5 

7.. £9 1 . £ 0 



RDI'J 

C 


MD GENERATORS 


TDP GEMERRTDR 

EDTH GEMEPRTDRS 

v .-yo 

CRL 

X -■■■' Y 3 

CRL 

X.--Y3 CRL 

£ . 5 0 

3 .97 

£ . 5 0 

4.13 

£.50 4.10 

3 '-i 

3 7 3 

4. 17 

3 . 43 

3.33 3.04 

4. 17 

£ . 15 

5. 0 0 

5. 05 

4.17 4. 0£ 

5. 0 0 

C. . ( c. 

5.33 

c ■ «’ c 

5.00 5.34 

CJ .-I ~j 

1 . 6y 

7. 71 

1 - 0 0 

5.83 3.37 

P.71 

0. 93 



7,71 1 . 36 



RDI..J 

D 


MD GENERATORS 


TDP GEMERRTDR 

EDTH GEMEPRTDRS 

X V 3 

CRL 

X--Y3 

CRL 

X Y 3 CRL 

1.67 

4. £4 

1 . 67 

3 . 96 

£ . i' l 3 • 3£ 

£. 71 

3. 48 

£.71 

3. 3£ 

3.54 3 . 4£ 

3. 54 

£ . 69 

3.54 

£.47 

4.37 7. 15 

4 . 37 

£. 4£ 

4 . 37 

3 . 5 1 

5 . £ 1 5.76 

5. £1 

3 . 1 5 



6 . £5 4. 0£ 

6. £5 

1 .71 

5. £1 

4.54 

7.5 0 0.33 

7.50 

0 . 9 0 

6. £5 

£ . 09 




7 . 5 0 

0. 93 
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of the location given by both 2-D prediction and the substitute gas, as was the case 
with 20° afterbody. The same merging of shocks, as with the 20° afterbody, is 
indicated by the data. 

As can be seen from all of these data, the heating rates in the combustion pro- 
ducts are very high. For the 20° afterbody shock free case (the configuration most 

closely resembling the actual flight article), calculated heating rates in the vicinity 

2 2 

of the nozzle exit are on the order of 300 cal/cm sec (53 BTU/ft sec) on a full 

2 

scale vehicle. Corresponding zero degree model rates (6000 cal/cm sec) are too 
high to measure. Near the intersection of the shocks from the cowl and side genera- s 
tors the rates go as high as 1358 cal/cm sec on the model, or 68 cal/cm sec on the cl e. 
full scale vehicle. 

NOZZLE SEPARATOR WAKES 

In order to study the effect of nozzle separators on the heat transfer distribution 
a series of experiments was made using the no. 2 substitute and combustion gases. 

The first of these experiments was made without the nozzle separators in order to 
extablish a basis for comparison. Heat transfer data were obtained at locations closer 
to the nozzle exit than in the previous work. These data, when plotted against Z/Y^, 

Fig. 79a, show a depression of heating rates for both gases toward the nozzle comers. 
We believe that these depressed rates indicate vorticies formed by nozzle comer 
flows and are probably responsible for the early departure from 2-D behavior noted 
in our = 8 experiments (Ref. 2). 

The flow on the internal nozzle side wall is generally at a higher pressure than 
the flow on the floor at the same streamwise station. When this occurs in a super- 
sonic flow, there must be a crossflow leakage in the streamwise corner boundary 
layer, from the sidewall out across the floor of the nozzle. The surface streamlines 
on the afterbody are therefore likely to be inclined toward the centerline at the span- 
wise corners of the nozzle exit. Since the Mach waves emanating from the ends of 
the nozzle side walls travel relative to the local streamlines, it seems likely that 
this process could produce the observed reduction in pressure below the two-dimen- 
sional case ahead of the Mach wave emanating from an undeflected streamline. This 
postulate is reinforced by the behavior of tests with extended side plates, in which 
the two-dimensional region was extended by the plates, but the infiltration of the 
lateral expansion was much more rapid at the end of the extended plates than it had 
been for the case of the nozzle alone. In this case the corner crossflow would get 
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Fig. 79 Heat Transfer Rate Distribution, 20° Afterbody, Q/Q(^ vs ' Z/Y 3 




much stronger along the plates because the floor (afterbody ramp) pressure is 
greatly reduced by the 20° turn at the nozzle exit. 

The data without separators also show a difference between the No. 2 substitute 
and combustion gas heating rates. The difference is apparent in Fig. 80 for the row 
of gauges at X/Y^ = 0. 781, a location just aft of the nozzle exit. For this first row 
of heat transfer instrumentation the combustion gas heating rates are lower than 
either the No. 1 substitute gas rates or the turbulent Price and Harris prediction. Our 
first thought on this phenomenon was that we had vibrational freezing of N^ with the 
combustion gas so we made calculations of that effect. Although the amount of energy 
in vibration is probably enough to account for the effect (30% of the gas enthalpy), the 
relaxation time for N 2 in the presence such large quantities of £^>0 is not long enough. 
The relaxation time would have to be on the order of 10 u sec rather than the calculated 
levels of 2 m sec). A second explanation of the observed data is that there is a 
difference in the tendency toward relaminarization of the turbulent boundary layer 
under the influence of the strong favorable pressure gradient at the intersection of 
the nozzle exit and the 20° afterbody. 

Several authors (Kearney, et al, Ref. 12) have found that the Stanton number 
of incompressible flows is reduced markedly by local zones of favorable pressure 
gradient. Kearney, et al, show that the reduction is truly local, in that for the case 
they studied the Stanton number returned rapidly to the flat plate level (as if there had 
been no upstream favorable gradient section) downstream of point at which the pressure 
gradient was again zero. This behavior is remarkably similar to that of the combus- 
tion gas relative to the substitute gas, the accepted critical value of the relaminariza- 

-v dp _0 

tion parameter (K =^3^ ~j^r) is 3. 7 x 10 . This was established by S. Kline in 

1967, since endorsed by many (e.g. , Ref. 12). At this critical value the "violent 
eruptions away from the wall near the edge of this viscous sublayer ceased when 
K = 3. 7 x 10” . " Even at values less than critical, skin friction and heat transfer 
levels have been shown to drop well below their corresponding flat plate levels (Refs. 

12, 14). 

Three new questions are posed by our measurements: 1) what is the effect of 
compressibility on relaminarization criteria? 2) what is the effect of an impulsive 
change in pressure level, rather than a finite section of finite pressure gradient? 

3) what is the effect of cooling or relaminarization? Without pretending to answer 
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ROW E NO GENERATOR 


L/1 




these questions, we do know that there is ample precedent for cooling of high Mach 
number boundary layers delaying transition, which is also a phenomenon related to 
the bursts of turbulence referred to above. We also know that the similitude re- 
quirements of our testing ensure that K is invariant among all three cases: flight, 
combustion and substitute gases, because K can be written as 

K = — d <pA0 

R d (X/D) 
e 

In view of the above discussion, it is reasonable to postulate that the substantial 
reduction in heating rate at the start of the afterbody expansion seen for the combus- 
tion gas is primarily due to the effect of wall cooling in enhancing the proven ability 
of favorable pressure gradients in reducing turbulent transport phenomena. It is 
probably due only in small part to slow relaxation of vibrational states, because 
the water vapor concentration is high enough to make that process more rapid than 
the observed decay. These postulates should be verified by future testing, as the 
correct understanding of cause and effect will have a large effect on the required 
thermal protection on a scramjet afterbody. 

The addition of nozzle separators has the effect of reducing the area ratio of 
the individual nozzle sections. This in turn increases the exit Mach number and lowers 
the exit density. All plotted separator heating rates have been scaled upward by the 
ratio of exit plane density without separators to exit plane density with separators so 
that comparisons can be made to the experiments without separators. The tabulated 
data have not been adjusted to account for this effect. It should be noted that the 
change in area ratio is different for the individual, nozzle channels. For instance, 
the center section of the nozzle with the separators placed to divide it into three 
sections has a density 18% lower than the outer sections. The data in Fig. 79b show 
the heating rate distribution versus spanwise distance (Z/Y^) for this configuration. 
Heating rates just downstream of the separators are increased about 20% over their 
value without separators while the center line values are decreased by nearly 20%. 

The depression of centerline rates disappears after X/Y x = 2.35, with heating rates 
returning to values at or above the rates with no segments. 

Fig. 81 a, b, c shows the heating rate distribution versus Z/Y 3 relative to 
nozzle separator for one nozzle separator. These data are scaled by the centerline 
heating rate at the respective values of X/Y 3 with no nozzle separators. They were 
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Fig. 81 Heat Transfer Rate Distribution, 20° Afterbody, Q/Qq^ vs Z/Y3 



taken by moving the separator to several slightly different spanwise locations, so that 
the fixed gauge locations could be at different displacements from the separator. For 
these tests the density at the exit of each nozzle section is different, causing the wake 
to turn in the -Z/Yg direction so it is not surprising that the heat transfer rates in 
the separator wake regions are not symmetric. The data at X/Y = 1. 76 and X/Y 

O *J 

= 2. 35 show a definite structure with a maximum heating rate between 20 and 40% 
higher than the centerline value. The data at X/Y^ = 0. 781 do not show a clear-cut 
structure, but are depressed below their levels withput separators. These observa- 
tions imply that the wake region is very narrow relative to the gauge width at upstream 
locations and broadens downstream to a resolvable width at the gauge spatial resolu- 
tion (see previous discussion under Instrumentation in section 2). 

The observed peaks of the heating distributions are not always at the separator 
location but further to the nozzle side. This shift may be caused by the wake turning 
or may be an artifact of the spatial resolution. The important point is that the mea- 
surements show a rapidly varying level of heating; a structure that is consistent with 
a strong streamwise vortex (or vortex-like crossflow) along the root corners of the 
separators, spreading slowly as it passes out over the afterbody surface. 

In order to investigate the effect of trailing edge bluntness one separator’s 
trailing edge was ground off so that it was 0. 508 mm wide, i. e. , 4 times the initial 
width. Figure 82 shows these data plotted against Z/Y relative to the separator 

O 

position and normalized as before. At X/Yg of 1. 56 the wake structure is slightly 
more pronounced than for the same location with the thin trailing edge, while at 
X/Yg = 2. 35 the rates are lower than in the previous case. The first observation 
is due to a wake that is wider relative to the gauge size. The lower rates at X/Y q = 
2.35 may be an indication that the maximum relative heating rate now occurs further 
upstream which would also be a function of increased wake size. The data at X/Y_ = 

O 

0.781 are also lower in this case and it not as logical to assume that this is caused 
by a blunt trailing edge. This effect is probably a relaminarization. If the boundary 
layer is on the verge of relaminarization any small stabilizing effect could cause it 
to relaminarize. 

The heat transfer data for these nozzle separator wake tests are tabulated as 
functions of X/Yg in Tables 21 through 24. 


78 ' 






TABLE 21 Q/QA SUBSTITUTE GAS 2 


EDTH SEGMENTS 


0. 78 

0.71 


1 . 57 

0 .56 


1 . 57 

0. 54 

RDM 

8. 74 

0. 45 


3.53 

0.41 


3. 53 

0 . 45 


4 . 3 1 

0. 57 


5 . 09 

0 . 35 


5.87 

0. 32 


EDTH SEGMENTS 



0. 78 

1 . 06 


1 . 57 

1.15 


er 

L_ • 

1 . 05 

RDM 

2 . 94 

0. 62 


"> • i"’ o 

0.41 


O . f' o 

0. 84 


4 . 5 0 

0 . 4:3 


4.50 

0. 43 


5. 29 

0 . 36 


6. 07 

0 . 39 


BOTH SEGMENTS 



0. 78 

0 . 6 6 


1 . 58 

0.64 


C* 

1 ■ V 1 

0. 62 

RDM 

3. 13 

0 . 92 


3.92 

0. 52 


4.70 

0 . 5 1 


5.48 

0 . 36 


6 . 26 

0 . 23 


7. 44 

0 . 2 1 


EDTH SEGMENTS 



0. 78 

0. 75 


0. 78 

0 . 96 


1 . 56 

0 . 9 1 

RDM 

O “» cr 

i— • 

0 . 84 


>! . • “! 3 

0. 54 


4.11 

0 . 69 


4. 1 1 

0.69 


5. 68 

0. 66 


6 . 46 

0 . 3 0 


EDTH SEGMENTS 



0. 78 

0.71 


0.78 

0. 75 


1 . 56 

0. 75 

RDM 

1 . 56 

U . ( 8 


-i er 

C. . 

0.59 


2 . 35 

0. 58 


3.52 

0.51 


4.31 

0. 45 


4 . 31 

0. 34 


5 . 8 7 

IJ . >! 6 


6 . 6 6 

0. 38 


6.66 

0. 34 



80 


TABLE 22 Q/QA COMBUSTION HEATING 




EDI. I 

fl 



HD SEGMENTS 


CEHTER 

SEGMENT 

BDTH SEGMENTS 


0 . 7y 

0. 63 

0. 78 

0 . 67 

0 . 78 

0.61 

1 . 57 

0 . 65 

1.57 

0. 62 

0. 78 

0 . 56 

1 . 57 

0. 65 

1 . 57 

0.62 

1 . 57 

0. 44 



ROW 

B 



HD SEGMENTS 


CEHTER 

SEGMENT 

BDTH SEGMENTS 


1 . 57 

1 . 03 

0 . 78 

!J . A b 

0. 78 

0 .67 

■“ ”| c* 

CL ■ 

0. 74 

1 . 57 

0. 77 

1 . 57 

1 . 09 

2.35 

0. 74 

2.35 

0. 79 

■-« cr 

a . 

0 . 6 1 



ROW 

C 



HO SEGMENTS 


CEHTER 

SEGMENT 

BDTH SEGMENTS 


0 . 78 

0. 85 

0. 78 

0 . 64 

0. 78 

0 . 49 

0. 78 

0. 79 

1 . 57 

1 . 48 

0.78 

0. 54 

1.57 

0. 99 

2 . 35 

0. 67 

1.57 

0. 52 

3. 13 

0. 83 

3 . 92 

0. 69 

1 . 57 

0.51 

4.70 

0. 62 

4.70 

0.71 

2.35 

0. 77 

4.70 

0. 57 

6 . 2 6 

0. 36 

CL m o 5 

0. 77 

5.48 

0. 52 

6 . 26 

0. 36 

3 . 1 3 

0 . 6 6 

6- . 26 

0 . 3 0 

6 . 26 

0. 36 

3.92 

0 . 6 0 



EDI.I 

D 



HD SEGMEHTS 


CEHTER 

SEGMENT 

BDTH SEGMEHTS 


0 . 7 8 

0. 69 

0 . 78 

0. 72 

0. 78 

0. 79 

1 . 56 

1.01 

1 . 56 

0 . 7 0 

1 . 56 

0. 75 

1 . 56 

1 . 04 

■j “i c r 

1.17 

2. 35 

1 . 06 

-”r cr 

d m *1' -I 

0. 73 

O . . O 

0 . 4 0 

•T; “• *”i 

0. 45 

•j ~J *=j 

0. 83 


0 . 4 0 


0. 45 



RDW 

E 



HD SEGMEHTS 


CEHTER 

SEGMENT 

BDTH SEGMENTS 


0. 78 

0 . 78 

0. 78 

0 . 65 

0. 78 

0 . 8 1 

0. 78 

0. 75 

1 . 56 

0. 67 

0. 78 

0. 77 

1 . 56 

0 . 9 1 

~i cr 

C. . 

0 . 3 9 

1 . 56 

0 . 8 1 

■“< O CT 

L_ m -~ T ■_* 

0 . 58 

2 . 3 5 

0. 39 

2.35 

0. 43 


81 



TABLE 23 Q/QA COMBUSTION GAS 





RDM 

ft 

:EGMENT 

9- . 

274 

SEGMENT 

9 — . 552 

0. 78 


0. 74 

0 . 78 

0 . 66 

1 . 57 


i j . 9 1 

0.78 

0 . 7 3 




ROM 

E 

:EGMEMT 

9-. 

2 i r' 4 

SEGMENT 

9 -.552 

0 . 78 


0. 68 

0. 78 

0. 62 

0. 78 


0. 58 

2.35 

0. 53 

2.35 


0 . 59 

2.35 

0 . 5 3 




ROM 

C 

EGMENT 

9-. 

8 r"' 4 

SEGMENT 

9 — . 55 2 

0 . 78 


0. 55 

0 . 78 

0. 55 

0. 78 


0. 53 

0. 78 

0. 54 

1 . 57 


0 . 76 

1.57 

0. 77 

2.35 


0. 67 

4.70 

0 . 63 

3.92 


0 . 56 

4.70 

0 . 6 1 

4 . r" 1.1 


0. 62 

5.48 

0. 34 

5.48 


0.41 

6 . 26 

0 . 3 0 

6 . 2 6 


0. 28 

6 . 26 

0 . 3 0 




ROM 

D 

EGMENT 

9 — . 

274 

SEGMENT 

9 -.552 

0. 78 


j j ,64 

0. 78 

i.i .65 

0. 78 


0. 63 

1 . 56 

0. 79 

1 . 57 


0. 87 

2.35 

0. 94 

1 . 57 


0. 96 


0 . 63 

2.35 


0 . 69 


0. 53 

-~i ■ j •“ 


0.55 

-j 

0. 53 



0. 54 

S 

0. 53 




ROM 

E 

EGMEMT 


274 

SEGMENT 

•~i c r ■-> 

■v . a. 

0. 78 


0 . 6 7 

0. 78 

0. 65 

0. 78 


0. 64 

1 . 56 

0 . 76 

1 . 56 


0 . 83 

1 . 56 

0 . 8 1 

1 . 56 


0 , 65 

1 . 56 

0 . 8 1 

2.35 


0. 37 

1 . 56 

0 . 8 1 


82 



TABLE 24 HEAT TRANSFER RATE SUBSTITUTE GAS 2 


r 




ROW 

FI- 

— 

- 






□ SEGMENT 


SEGMENT 

S' 


274 

SEGMENT 

S' 

-. 552 


0.78 

0. 79 

0. 73 


0 , 

73 

0. 

i y 



ij . 6 6 

0. 73 

0 . 32 

2.74 


0. 

55 

2. 

74 



0 . 5 0 

1.57 

U ■ r »■’ 

3.53 


ij . 

63 

“i 

53 



0 . 6 0 

1.57 

U . >■' cl 

4.31 


o. 

43 

4. 

31 



0 . 37 

2.74 

0. 62 

5. 09 


o. 

9 7 

5. 

09 



0 . 33 

3. 53 

0 . 73 

5.37 


o . 

6 

5. 

37 



0 . 32 

C •~I 

. “i m .ji -Z' 

0.75 

5. 37 


0. 

36 

5. 

z« r' 



0. 32 

4.31 

0. 47 

5.37 


l”j 

36 

5 ■ 

37 



0. 32 

4.31 

0. 59 

5.37 


0. 

36 

5. 

"i "7 
J l 



0. 32 

5. 09 

0 . 49 

5.37 



36* 

5. 

z« r 



0. 32 

3 . 66 

0 . 4 1 

5.37 


0. 

36 

5. 

Z« f 



0. 32 

7. 44 

0. 32 

5.37 


I-I 

36 

5. 

-•7 



0 . 32 

.“i •— i 

m C.C. 

0 . 25 

5.37 


0 . 

36 

5. 

37 



0 . 32 

!9 . 0 1 

0 . 2 1 

5.37 


0 . 

36 

5. 

Z* i 



0 . 32 



ROW 

B- 

— 

- 






□ SEGMENT 


SEGMENT 

S< 


274 

SEGMENT 

S' 

-.552 


0.73 

0 . 93 

0 . 73 



6 3 

0 . 

r o 



0 . 69 

1.57 

0 . 9 1 

1 . 57 


n . 

79 

1 . 

57 



0 . 39 

2. 35 

0 . 3 1 

-* ‘ “I ^ 



31 


35 



0. 87 

2.94 

0.71 

2 . 94 


0 . 

i** O 

|” m 

34 



0 . 69 

•Z* m i a 

0. 74 

4.50 


0 . 

39 


78 



0 . 7 0 

4.50 

0. 32 

5-. 29 



47 

4. 

5 0 



0. 38 

4.50 

0 . 34 

6 . 07 


o . 

O *“* 

5. 

29 



0 . 5 0 

5. 29 

0 . 5 0 

6 . 0 7 


0 . 

O •“ 

6 . 

07 



0. 38 



ROW 

C- 

— 

- 






ID SEGMENT 


SEGMENT 

3 


2 r" 4 

SEGMENT 

S' 

-.552 


0 . 73 

1 . 15 

2.35 


0 . 

91 

0 . 

r 'Z' 



0 . 74 

1 . 56 

1 . 02 

3 . 1 3 



99 

1 . 

56 



1 . 0 0 

2.35 

0 . 34 

3. 92 



7 0 

C- • 

-,CT 



0 . 95 

3 . 1 3 

0 . 9 1 

4.7 0 


0 . 

7 0 

o 

32 



0. 69 

• - 1 . 9 c! 

0 . 57 

5 . 43 



51 

4. 

7 0 



0. 72 

4.70 

0 . 5 0 

6 . 26 



“i 

5. 

48 



0. 52 

5.43 

0. 29 

7.44 



31 

6 . 

26 



0. 35 

6 . 26 

0. 23 

7.44 



31 

!'" ■ 

44 



0 . 3 0 



ROW 

II- 

— 

- 






□ SEGMENT 


SEGMENT 

S' 


8 >■' 4 

SEGMENT 

S' 

ETC O 
m ■— 1 ■-.* L— 


0 . 73 

1 . 09 

1 . 56 


1 . 

24 

0 . 

7 ;Z; 



0. 42 

1.56 

0 . 96 

in’ ■ -z*5 


1. 

13 

1 . 

56 



1 . 04 

,L 

■£* m -Z» 

0. 89 

O . >!>! 


0 . 

44 

c ■ 

Z' J 



1 . 09 

4. 1 1 

0 . 62 

4. 1 1 


0. 

65 

3. 

z'-z* 



0. 86 

4. 39 

0 . 44 

4 . 39 


0 . 

47 

4. 

1 1 



0 . 68 

5.63 

0. 43 

5 . 63 


0. 

72 

4 . 

39 



0 . 52 

6. 46 

0 . 23 

6.46 


o. 

45 

5 . 

Sy 



0. 62 

7.44 

0 . 3 1 

6. 46 


0 . 

45 

5. 

S3 



0. 62 



ROW 

E- 

— 

- 






0 SEGMENT 


SEGMENT 

3 


8 r' 4 

SEGMENT 

S' 

- . 552 


j . 56 

0 . 6 7 

|“| "p 1 0 


0. 

y 6 

0. 

r'* O 



0. 37 

It! 

0 . 66 

1 . 56 


0. 

86 

1. 

56 



0. S3 

4.31 

0. 51 

2. 35 


0. 

74 

C- m 

35 



0. 63 

5. 09 

0 . 38 

3. 52 


0 . 

31 

o 

52 



0, 33 

6 . 6 6 

0 . 34 

4.31 


0. 

49 

4. 

31 



0 . 52 

7.44 

0. 34 

5. 09 


0. 

43 

5. 

09 



0.43 

7. 44 

0 . 34 

5 . 37 


0. 

39 

5. 

37 



0. 32 

7.44 

0 . 34 

6 . 66< 


0. 

51 

6. 

66 



0. 43 


83 





COWL SHOCK IMPINGEMENT LOCATION 


The calculation of the cowl shock impingement location showed very little change 
with gas composition. The data in Fig. 83 show the heat transfer distributions for 
both substitute gases with the cowl shock generator. These data agree with the calcu- 
lations and show essentially the same shock location. It should be remembered that 
the No. 2 substitute gas experiments were made with an extended nozzle side to in- 
crease the 2-D region and that this accounts for some of the apparent scatter. 

As a result of those experiments and the calculations it is apparent that the aft 
cowl shock location for the combustion gas is caused either by a relaxation effect 
(about lyu sec of relaxation time required) or a shock boundary layer interaction. 
Which of these two candidates is correct could be determined by changing the combus- 
tion gas temperature to see if the cowl shock impingement location moves. 

Explaining the difference in shock location between combustion and substitute 
gases is important, even though it is small in magnitude. We must understand the 
limitations of the substitute gas method if it is to be used with confidence in develop- 
ment testing. Since it is clear from the above that the differences are not due to 
small differences in substitute gas properties, we look to boundary layer and non- 
equilibrium effects for their origin. We believe the answer lies in the effect of wall 
cooling on shock boundary layer interaction, but it is also possible that there is an 
effect of finite vibrational relaxation time in the combustion gas that is not present 
in the substitute gas. The impingement of an incoming shock on a boundary layer 
has a well characterized pressure distribution, the streamwise extent of which de- 
pends on the boundary layer and shock characteristics. The difference in wall cooling 
between substitute and combustion gas would act to make a separated zone around the 
impingement more compact in the combustion case by increasing the density in the 
separated zone. The pressure and heat transfer rises would therefore not extend 
as far upstream as in the substitute gas case. This explanation is plausible because 
the shifts in position of the heat transfer rise between combustion and substitute gases 
are small compared to the streamwise distance over which the rise takes place. Al- 
though the cooling effect on separated flows has been observed in previous investiga- 
tions of cooled separating flows at high Mach number (e.g. , Kaufman, et al, Ref. 13), 
it is quite speculative in this case. Kaufman, et al, found that "In general, cooling 
the wall delayed separation and reduced its extent, except when the cooling was 


84 



VD/D 






t 


i 




t* 


X/Y3 


Fig. 83 Heat Transfer Rate Distribution, 20° Afterbody, with Top (Cowl) Generators, All Rows 



sufficient to change the character of the boundary layer. ” Their tests were for a con- 
figuration very similar to the present shock generator, so the probability also exists 
that some of the shift in shock location is attributable to the shock generator boundary 
layer interaction. 


6. CONCLUSIONS AND RECOMMENDATIONS 

The detonation tube has provided a good simulation of the exhaust flow from a 
hydrogen/ oxygen scramjet engine at a free stream Mach number of 6. Pressure 
distributions obtained on a cowl and afterbody model with the flow of simulated com- 
bustion products and the flow from a substitute gas mixture of 50 percent Argon and 
50 percent Freon 13B1 were in good agreement in the two-dimensional regions of the 
flow. 

Extensive heat transfer measurements were in very good agreement with theo- 
retical calculations for a turbulent boundary layer. This is in contrast to our measure- 
ments at the M qq - 8 flight condition (Ref. 2), where the measurements compared well 
with laminar boundary layer calculations. The extremely high heating rates mea- 
sured in the present tests should be well noted. Because these data were obtained at 
full scale Reynolds number and total enthalpy, they can be directly converted to full 
scale vehicle values by dividing by 20, the model scale factor. They show full scale 
heating rates the order of 18 eal/cm sec at a distance of 3 combustor exit heights 

downstream of the combustor exit on the 20° afterbody (as close to the exit as we 

2 

measured) and close to 70 cal/ cm sec further downstream where two artificially 
generated shocks crossed. Because the actual scramjet vehicle will have several 
engine modules and the module separators will generate wakes and crossing shocks, 
we recommended additional testing, more closely duplicating actual vehicle geometry, 
to determine the effects of such disturbances. Our additional testing showed that 
separator wakes do have significant effects on the afterbody heating distributions. 

There were local increases in heating rate as much as 60% above the centerline 
value with no separators. The peaks were strong gradients in the heating rate 
which must be accounted for in vehicle design. Furthermore, due to limited spatial 
resolution, the actual peak heating rates on gradients may be much higher than the 
indicated values. The experiments suggest that there are vortical crossflows induced 
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at the corners and that these are responsible for the large changes in the heating 
distribution. The combination of these vortices and the separator wakes make the 
actual vehicle afterbody heating distributions difficult to predict. We recommend 
additional testing with a model equipped with heat transfer instrumentation specially 
made to give high spatial resolution in both streamwise and spanwise directions. The 
model should be close to the actual vehicle geometry as possible. 

The experiments performed under the extended contract suggest that a greater 
tendency toward re laminar ization of the turbulent boundary layer occurs with rapid 
expansion of the combustion gas than with the substitute gas. The difference is 
probably caused by the stabilizing effect of the relatively colder wall in the combustion 
case. Additional experimentation and calculation along these lines could show the way 
to reduce the afterbody heating rates on the actual vehicle by confirming that this is 
the real cause of observed heating rate differences. 

One anomaly left unexplained by the first tests at M w = 6 was the slight difference 
in the cowl shock impingement location on the afterbody for the substitute gas com- 
pared to the combustion gas. This particular substitute gas mixture was chosen be- 
cause it was the same as that to be used in a concurrent test program at the NASA 
Langley Research Center, and our two-dimensional M-O-C calculations showed that 
it very closely matched the calculated flow field using equilibrium combustion gas 
products. We recommended that additional calculations and tests be made with slightly 
modified substitute gas mixtures in order to bring the shock positions into more per- 
fect agreement. 

The results of additional experiments and calculations showed that the cowl 
shock impingement location is only very weakly related to the substitute gas composi- 
tion. Therefore, the anomaly is not caused by differences in pressure ratio or shock 
angle. The aft arrival of the combustion gas cowl shock must be caused by either a 
relaxation effect in the combustion gas or shock boundary layer interaction, which 
can differ in the two cases because of the wall cooling effect. 

We encountered a problem in this test program when we tried to measure impact 
(Pitot) pressures in the exhaust flows. The purpose of the measurements was to cali- 
brate the flow to further verify the simulation technique. Basically, our probes were 
too large for the gradients and size of the flow fields and meaningful measurements 
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couldn’t be obtained. For diagnosing small flow fields such as these we would recom- 
mend that non-flow-disturbing instrumentation techniques, such as laser-Doppler- 
velocimetry, be used. 
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